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Abstract —When multiple radio-frequency sources are con¬ 
nected to multiple loads through a passive multiport matching 
network, perfect power transfer to the loads across all frequencies 
is generally impossible. In this two-part paper, we provide 
analyses of bandwidth over which power transfer is possible. 
Our principal tools include broadband multiport matching upper 
bounds, presented herein, on the integral over all frequency of 
the logarithm of a suitably defined power loss ratio. In general, 
the larger the integral, the larger the bandwidth over which 
power transfer can be accomplished. We apply these bounds in 
several ways: We show how the number of sources and loads, 
and the coupling between loads, affect achievable bandwidth. We 
analyze the bandwidth of networks constrained to have certain 
architectures. We characterize systems whose bandwidths scale 
as the ratio between the numbers of loads and sources. 

The first part of the paper presents the bounds and uses them 
to analyze loads whose frequency responses can be represented by 
analytical circuit models. The second part analyzes the bandwidth 
of realistic loads whose frequency responses are available numer¬ 
ically. We provide applications to wireless transmitters where the 
loads are antennas being driven by amplifiers. The derivations 
of the bounds are also Included. 

Index Terms —Bandwidth, Bode-Fano bounds, broadband 
matching bounds, non-reciprocal networks, passive matching 
networks, radio-frequency coupling 

I. Introduction 

Analyses of the bandwidth over which power can be 
transferred in radio-frequency (RF) systems are often limited 
to a single source and load because of the complexity in 
manipulating multiport matching networks and multiple loads. 
Factors that contribute to the complexity include defining 
appropriate measures of bandwidth when there are many 
sources and loads, and the difficulty of analyzing coupling 
between loads. We propose methods of analysis that utilize 
broadband performance bounds applicable to a wide class 
of passive networks and an arbitrary number of sources and 
dissipative loads. 

The ability to transfer power from sources to loads relies, 
in part, on the ability to match the impedance of the sources 
to the frequency-dependent impedance Zi(jw) of the loads 
over a broad frequency range. Bandwidth upper bounds are 
of great help in determining the best achievable bandwidth 
performance for a given load. Classical Bode-Fano results [1], 
[2] on the integral of the logarithm of the reflection coefficient 
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can be used for such bounds when there is a single source 
and load. When there are multiple loads, analyses are often 
limited to special cases. For example, bandwidth bounds for 
loads with some specific structures are discussed in [3]-[6], 
and examples of analyses of multiport systems include [7], [8]. 
Often, multiple reflection coefficients are defined and analyzed 
separately using scalar Bode-Fano theory. However, as shown 
in [9], a physically-meaningful single reflection coefficient can 
be defined and analyzed when there are N arbitrarily coupled 
loads driven by N sources {N > 1). 

In this two-part paper, we present a bandwidth analysis 
of matched multiport RF systems that builds on bandwidth 
bounds in [9]. The first part presents the bounds and applies 
them to systems that can be expressed in closed form. The 
second part provides proofs of the bounds and applies them 
to systems whose scattering parameters are expressed numer¬ 
ically. The bounds in [9] apply to loads that are modeled as 
perfect reflectors as w —> oo. We extend these results, and 
present bounds that apply to loads that are reflectors at any 
frequency, including w = 0. We allow the matching network 
and loads to be non-reciprocal. The network can also be lossy. 
We permit the number of sources and loads to be unequal. 

By applying bandwidth bounds, we demonstrate how the 
number of sources, loads, and the coupling between loads 
affect the achievable bandwidth of a matched multiport system. 
We prove that bandwidth bounds generally scale as N/M, 
where M is the number of sources and N is the number of 
loads. This result also holds in the presence of coupling, as 
long as it is not “too strong”. This suggests that unlimited 
bandwidth is theoretically achievable by simply adding more 
loads for a fixed number of sources. As is shown, both the 
loads and the network architecture play an important role in 
achieving linear-in-7V performance of the overall system, for 
a given M. 

We also propose a bandwidth analysis for situations where a 
portion of the network is constrained to have a certain structure 
while other portions are unconstrained. This situation occurs in 
beamforming applications since a beamforming antenna array 
can be thought of as N loads driven with prescribed amplitude 
and phase relationships by a single source. 

The basic premise of broadband matching is that when 
a source and load are connected to each other, even if the 
reflection coefficient is made small at a design frequency 
w = ujd, it is generally not small for all oj [10]. When there 
are multiple sources and loads, there is no single reflection 
coefficient since power sent from source i may, through 
coupling, return to source j ^ i. In [9], a definition of a 
multiport reflection coefficient that takes this phenomenon 
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into account is used to derive bounds on the ability to match 
multiple sources and loads over all w with a lossless network. 
We expand this definition to include lossy networks. 

Of particular interest is the application to loads that are 
closely-spaced antennas, such as may be found in multiple- 
input multiple-output (MIMO) communication systems. The 
“densification” of portable wireless communication devices, 
including cellular telephones, with multiple transmitter and 
receiver chains in close proximity, makes coupling difficult 
to avoid. Furthermore, there are situations where there are 
more antennas than RF chains [11]. Our analysis methods 
quantify the bandwidth attainable in the transmitters of these 
MIMO systems, where the RF amplifiers are treated as sources 
and the coupled antennas are the loads. Part II, in particular, 
shows how realistic antennas are modeled to obtain accurate 
bandwidth results. 

We consider an RF system where M sources drive N 
loads through an arbitrary passive {M -b A^)-port matching 
network. The M input ports on the multiport network connect 
to the sources, and N output ports connect to the loads. 
No relationship between M and N is assumed. Our goal 
is to transfer maximum power from sources with known 
characteristic impedances to loads with known frequency- 
dependent impedances. The loads are dissipative and poten¬ 
tially non-reciprocal; the network can be lossy and also non¬ 
reciprocal. The quality of match between the sources and the 
loads at a frequency w is then determined by the power lost 
either because it is returned to the sources or because it is 
dissipated in the network. We derive and utilize bounds on 
this quality metric when an arbitrary passive (M -b A^)-port 
network is used. Our bounds can readily be calculated from 
the frequency-dependent S-matrix of the loads. 

We do not consider noise in our analysis, and hence our 
methods apply more to transmitters driving loads than to 
receivers, where the noise figures of the input amplifiers play 
an important role. We are interested in the “input bandwidth” 
of the loads and do not evaluate the ability of the loads to use 
their input power efficiently. 

When the loads are coupled, there is only a nominal 
association of source i with load i for i = 1,...,7V since 
source i potentially also stimulates load j ^ i \f the two 
loads are coupled. In minimizing the amount of lost power, 
a matching network between the sources and loads therefore 
could connect any source with any of the loads, and converse¬ 
ly; examples include decoupling networks [12]-[15], which 
ensure all power from the sources is delivered to the loads at 
a design frequency. An effective network prevents reflection by 
decoupling the loads from each other over as wide a frequency 
range as possible. Generally, for a given design frequency 
uJd, there are frequencies oji < ujd and > uid, where 
the fractional power delivered to the loads falls below some 
prescribed threshold. The larger 022 — wi is, the larger the 
bandwidth. 

Techniques involving active non-Foster circuits [16], [17] 
and adaptive matching [18]-[21] are not discussed. Non-Foster 
circuits such as those used to realize negative capacitance 
and inductance values can theoretically cancel the reactive 
Foster behavior in the loads, and achieve an impedance match 



Fig. 1. An RF system where M uncorrelated sources drive N loads having S- 
matrix Sl through a passive (M+A^)-port matching network S (the complex- 
frequency argument s is omitted). We use Slm to denote the S-matrix as 
seen at the input of the network, and Sq = S 22 to denote the S-matrix as 
seen at the output. 


between the sources and the loads over a wider frequency 
band than passive networks. Adaptive matching uses tunable 
capacitors and inductors to match antennas whose impedance 
is affected by the environment or carrier frequency. 

The basic model assumptions are given in Section II. The 
definitions and notations we use are similar to those in [9]. The 
principal bounds we use are presented in Section III. Details 
on the mechanics of how to apply these bounds are provided 
within the subsections of Section III. Analytical applications 
of the results appear in Section IV, where scaling laws are 
derived and network architectures are analyzed. Section VI 
concludes. 

The results presented herein assume the loads can be accu¬ 
rately modeled by rational functions of frequency. Since not 
all loads are necessarily rational. Part II examines how rational 
functions can be used as approximations for arbitrary loads. 
We present examples of how to use the bounds in practice. 
Proofs of all the results are also given. 

II. Problem Definition AND Notations 
A. System description 

Figure 1 shows an RF system where N dissipative loads 
are driven by M decoupled sources with real impedance Zq, 
the characteristic impedance of the system. Let Sl{s) be the 
N X N S-matrix of the dissipative loads as a function of the 
complex frequency s = a+juj, where tr and w are real. Sl{s) 
is obtained by extending the S-matrix of the loads as 

a function of the angular frequency w to the whole complex 
plane (WCP). Mathematically, Sl{s) can be thought of as the 
transfer function between the N x 1 incident and reflected 
waves 026 ®* and 626 ®*, where t represents time. Therefore, 
we have 62 ( 5 ) = S'L(s)a 2 (s). If the loads are reciprocal then 
Sl{s) is symmetric. If the loads are coupled, at least one off- 
diagonal entry of Sl{s) is non-zero. The impedance matrix 
ZlU^) of the loads can be obtained by ZL^joj) = Zq{I + 
SL{juj)){I - SlUuj))-^. 

The M sources and N loads are matched by inserting a 
passive (M + V)-port network between them, as indicated in 
Figure 1. The M input ports of the network are connected 
to the sources, and the N output ports are connected to the 
loads. The network is not necessarily lossless or reciprocal, so 
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we allow standard passive capacitive and inductive elements as 
well as non-reciprocal ferromagnetic components in its design. 

Let S{s) be the (M -|- N) x (M -|- N) S-matrix of the 
multiport network as a function of s, partitioned as in Figure 1 

M N 

N \S2l{s) S22{S)J 

Let Slm{s) denote the M x M S-matrix of the cascade of 
the network and the loads, and Sa{s) be the N x N S-matrix 
seen from the output ports of the network. Since the input 
to the network is terminated by sources with characteristic 
impedance Zq, it follows that 5'g(s) = S 22 {s), and 

SLM{s) = Sn{s) + Si 2 is)SLis){I - Sg{s)Sl{s))-^S 2 i{s). 

( 2 ) 

As in Figure 1, let M x 1 vector ai(s) be the incident 
wave to the network. Then the reflected wave is &i(s) = 
SLM{s)ai{s). Since the network is potentially lossy, the 
incident power may be reflected to the sources (return loss), 
or dissipated in the network (insertion loss). We now define a 
measure of performance that includes both of these effects. 

B. Power loss ratio 

Our measure of performance of a matching network be¬ 
tween the sources and loads is given by the fraction of source 
power that is not transferred to the loads, as a function of 
frequency. At frequency s = joj, the total instantaneous power 
from the M sources is ||ai(ja;)|p, and the total instantaneous 
power delivered to the N loads is ||a 2 (jw)|p —|| 62 (jw)|p > 0. 
The power lost due to dissipation and reflection is therefore 
||ai(j.;)f-(||a 2 (j..)f-|! 62 (ju;)f). 

Definition 1: The power loss ratio at frequency joj is the 
ratio between the expected power loss and the expected total 
incident power at jw: 

2 , , _ E[||ai(j.3)f - {\\fi2{3w)r - Whiju^W)] ... 

E||ai(ja;)P ’ 

where the expectation is over the random input signals. 

By convention, when we use r{uj) we mean the positive 
square root of (3), and by construction, 0 < r{uj) < 1 where 
values close to zero indicate that little source power is being 
lost and therefore most of it is being delivered to the loads. 
Values close to one indicate most of the power is lost to 
dissipation or reflection. We note that r(w) = 0 means that the 
loads are perfectly matched and decoupled from one another. 
When the matching network S{s) is lossless, the power loss 
ratio is equivalent to the power reflection ratio defined in [9]. 
By considering power dissipation in the network in addition 
to power reflection by the network, we handle situations that 
appear favorable because the reflected energy is low, but are 
actually unfavorable because the network (rather than the load) 
is dissipating the incident power. This issue is raised in [22]. 


C. Experimental measurement of r(uj) 

Because (3) includes an expectation, we devote a few words 
to the measurement of r{uj) in practice. The incident signals 
from the sources are assumed to have independent random 
phases and instantaneous powers at each frequency, and the 
expectation is taken over all such uncorrelated incident signals. 
Nevertheless, the expectation is unnecessary when M = 1 
since the amplitude and phase of a single source does not affect 
the power loss ratio. In this case, 1 — r^(w) is equivalent to the 
standard transducer power gain [8], [23], and the expectations 
in the numerator and denominator of (3) may be dropped. 

When M > 1 the expectations play the important role 
of averaging over all amplitude and phase combinations of 
the input signals. Its importance can be demonstrated by 
examining M = 2, where it is well known that even-mode 
(in-phase) and odd-mode (out-of-phase) signals can elicit very 
different reflective responses from a two-port system. Since the 
sources are uncorrelated, no preference is given to even or odd 
mode signals, and the expectation eliminates mode dependence 
by averaging over both of them. Hence, r(a;) can be thought of 
as an “average loss” experienced when the loads are stimulated 
by uncorrelated sources. 

The value of r^(w) can be experimentally measured by 
forming the sample average of the numerator |jai(ju;)|p — 
(||d 2 (ja;)|p — \\b 2 {juj)\f) for a variety of inputs di(ja;), and 
taking the ratio of this average to the sample average of the 
denominator |jdi(jw)|p. This ratio of averages converges to 
r^(w) as more measurements are taken with all possible source 
phase and amplitude combinations. 

D. Definition of bandwidth 

The network S{s) should be constructed to make r{uj) as 
small as possible over a prescribed bandwidth, or make the 
bandwidth as wide as possible for a prescribed threshold. 
Usually, bandwidth is measured in the vicinity of a design 
frequency, which we denote as A decoupling network [15] 
enforces r{uj) = 0 at a; = ujd- We can define the bandwidth 
of the combined network and loads using (3). 

Definition 2: The bandwidth is the frequency range for 
which r{uj) is no greater than a threshold r > 0 in the vicinity 
of a design frequency ujd'. 

u;Bw(T,Wd) = max a; 2 -a;i. (4) 

a;i<CL>d<a;2 

[cji , 102 ] 

Let the elements of di{juj), representing the incident signal 
from M decoupled sources at frequency juj, have equal ex¬ 
pected square-magnitude over all frequency, and have uniform¬ 
ly distributed random phases in [0, 27r) that are independent 
of the amplitudes and each other. Then (3) yields 

2 / ' 1 ^ 1 _ Etrjaf (ja;)(/ - gf {juj)SL{ju:))d 2 {ju:)} 

Etr{af (ja;)ai(jw)} 

Because the phases are independent and uniformly distributed, 
E[ai{juj)d^ {juj)] is a multiple of the identity matrix. We 
apply a 2 (jw) = (/ - S'G(jw)5'L(jw))"iS'2i(jw)ai(jw) to 
obtain 

r2(w) = 1- 



4 


- SgSl)-^{I - Sl){I - SgSl)-^S2i} 

M ’ ^ ’ 

where tr( ) denotes trace, ^ denotes Hermitian transpose; the 
frequency argument joj is omitted on the right-hand side of 

(5). 

When M > N, > 1 — N/M since the matrix 

inside the trace on the right-hand side of (5) is a rank-A^ 
positive semidefinite matrix whose eigenvalues are less than 
or equal to one. Hence, a certain fraction of the source power 
is always reflected or absorbed, and letting r < 
in (4) always obtains zero bandwidth. We therefore assume 
- N/M < T < 1 for M > iV; for M < TV, we have 
0 < r < 1. 

Note that “ganging” amplifiers through couplers to attain 
high output power with a single load appears to be an example 
where M > N. But such ganged sources are correlated 
since each carries the same signal, possibly differing only 
in a constant relative phase or amplitude. Since we assume 
uncorrelated sources, ganged amplifiers (and other correlated 
sources) should be considered as a single source to apply our 
results. 

Although the incident signals from the sources are assumed 
to be uncorrelated, the output of the matching network will 
have correlated components when M < N since any network 
driving all N loads will necessarily derive its signals from the 
M sources. As an example, a beamformer can be modeled 
as M = 1 source driving > 1 loads in a fixed phase 
relationship. We consider beamforming in Section V-A. 

E. Properties of S-matrices 

We briefly summarize some properties of S-matrices since 
they play an important role. Passive real networks have S- 
matrices that are real-rational, Hurwitzian, and bounded; this 
applies, for example, to Sl{s) and S{s). The definitions 
of these terms can be found in [23]-[25]. We also employ 
definitions of poles and zeros of rational matrices that are 
widely used in multivariable control theory [26]. For an 
arbitrary rational matrix A(s): 

• Poles: are the roots of the pole polynomial of A(s), where 
this polynomial is the monic least common multiple of 
the denominators of all minors of all dimensions of A(s). 

• Zeros: are the roots of the zero polynomial of A(s), 
where this polynomial is the monic greatest common 
divisor of the numerators of all minors of dimension L, 
and L is the normal rank of A(s). It is assumed that these 
minors have the pole polynomial as their denominators. 

The normal rank of a matrix is its maximum rank among 
all s G C; the use of “normal” refers to the rank almost 
everywhere in the complex plane [24, A.70]. We use LHP 
to denote the left-half complex plane (Re{s} < 0) and RHP 
to denote the right-half (Re{s} > 0). Throughout, we use ^ 
and ZL,i, i = 1,2,... to represent the poles and zeros over 
the WCP (whole complex plane) of Sl{s). Since Sl{s) is 
Hurwitzian, it has no RHP poles. 

We assume I — S]^{—s)Sl{s) is full normal rank. This is 
satisfied if the loads are strictly dissipative, so that there is no 
combination of load signals that is completely reflected for all 


s. On the other hand, we also assume that there exists an sq 
with Re{so} > 0 such that 

SI{-so)Sl{so) = I. (6) 

In general, sq is arbitrary and can be infinite. 

When So is purely imaginary sq = jujQ, (6) has the 
interpretation of modeling the loads as perfect reflectors at 
frequency ojq since S'£{—joj) = S^{juj) and therefore the 
singular values of Sj^ijoJo) are all unity. Because Sg{s) is a 
bounded matrix, (6) is equivalent to | det5'i(ju;o)| = 1. 

When So has positive real part, the physical interpretation 
of (6) is elusive but we still use the terminology “perfect 
reflector” at so. Examples of load structures with such an sq 
are given in Section IV-B2. 

For any so, (6) must also hold if we replace sq with 
—So. Therefore, without loss of generality, we only consider 
Re{so} > 0. If there are multiple distinct values of sq for 
which (6) holds then the bounds presented herein apply to each 
value separately. We therefore consider only a single distinct 
So¬ 


ul Broadband Matching Bounds 

We present the principal bounds used in this paper. They are 
derived in detail in Part II, but knowledge of the derivations 
is not needed to apply the bounds. A description of how to 
use these bounds appears in Section III-B. Conditions for 
achieving equality are presented in Section III-C. The bounds 
depend on zeros and poles of S'l(s), and techniques to obtain 
these are discussed in Section III-D and Part IT 


A. Principal bounds 

Bound 1: For sq = jwo, 

(coo - Lo)-^ + {loo + u;)-^ 

Jo 2 


log 


r(w) 


duj 


^ +'^{ZL,^+julo) ^ 


(7) 


Bound 2: For Re{so} > 0, 

Re[(so - -f (so + 


log- 


1 


-'0 


detS'L(so) 


2 r{u) 

nj(so+ ZL,i) 


doj 


n*(so - ZL,i) 


Bound 3: For sq = oo, 

/■“, 1 , 

/ log ——aw < 

Jo ^ ^(‘^) 2M 




zla 


( 8 ) 

(9) 


Bound 1 is useful for loads that are modeled as electrically 
open or short at some frequency jujQ. For example, some 
antennas are capacitive relative to ground when they are 
“electrically small” compared to the signal wavelength. Thus, 
they are effectively an open circuit at sq = 0; equivalently 
Sl{0)=I. 

Bound 2 applies to loads that are modeled as a mixture of 
resistive and reactive components [27]. An example of this is 
given in Section IV-B2. 
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TABLE I 

Forms of f{jio) and B in (10) for different locations of so, 

WHERE ARE THE POLES AND ZEROS OF Sl{s). 


so = foJo 

/(co) = |[(aJo - co) 2 + 2] 

B [Ei(PL,i + Ei(2L,i 

Re{so} > 0 

/(to) = iRe[(so -juj) 1 + (so +juj) 1] 
B= -^log|detS^(so)-ng^| 


SO = 00 

f{‘o) = 1 

i PL,i + Xli ^L,i) 



Bounds 2 and 3 can be applied in a similar fashion. If 
S'f (oo)S'l(oo) = /, (9) gives 


UJ2 — UJi < 


log(l/r)’ 


( 12 ) 


where is the right-hand side of (9). This gives us a direct 
bound on the bandwidth wbw(t,W d) achievable for all Wd. 
Equations (11) and (12) are complementary in that both can 
be in force simultaneously. 

We now discuss conditions under which these bounds can 
be achieved. 


Bound 3 applies to loads that are modeled as open or short 
circuits at infinite frequency. The classical model of a parallel 
resistive and capacitive load used to demonstrate the Bode- 
Fano bound falls into this category. A bound similar to (9) is 
presented in [9]; however, the bound in [9] requires M = N 
and the matching network to be lossless. 

In all three cases, the bounds have the form 

fi^^) log -^duj < B, (10) 

r(a;) 

where /(w) > 0. The form of /(w) depends on the location 
of So, and the computation of B depends also on the poles 
and zeros of Sl{s). Because 0 < r{uj) < 1, we have 
log(l/r(w)) > 0. Hence, /(w) log(l/r(u;)) > 0 for any oj. 
Clearly, B must be positive as well. We use (10) to generically 
indicate any of (7)-(9). The number of sources appears as M 
in the denominators of all three bounds. The forms of /(w) 
and B are summarized in Table I. 



B. How to use bounds 


Suppose that we wish to assess the achievable bandwidth 
[uji,uj 2 ] of a set of loads, where uji < UJ 2 - Our measure of 
achievability is that for some threshold t > 0, the overall 
system should obey r{uj) < r for w € [ 001 , 002 ]- Hence the 
combined multiport network and loads reflects (or absorbs) 
no more than t within the passband. We assume that S'l(s) 
is available to us (we have more to say about this in Section 
III-D) and we would like to know wbw(t, Wd) defined in (4). 

Suppose the loads obey S'J(0 )S'l(0) = / so that Bound 1 
(first row of Table I) with wq = 0 applies to any passive 
network used for these loads. Let the right-hand side of (7) be 
denoted Bi > 0, which depends only on S'l(s). Then 


Bi > 


j ^ log } duo > log - 
r{oo) T 


= log- 

T 



The first inequality applies to any network, while the second 
inequality applies to a network with the desired passband 
characteristics. Hence, 


1 I ^ Bi 

LOi U 02 ~ log(l/r)' 


( 11 ) 


Clearly, this inequality imposes a constraint on (wi, 002 ) pairs. 


C. Conditions for equality 

We distinguish between conditions for equality in (7)-(9) 
and conditions for equality in (11), (12). The former apply 
to any passive real network and the latter apply to networks 
with particular passband characteristics. Bounds (7)-(9) have 
a common set of conditions for achieving equality: 

1) 52r(s)^Ji(-s) + Sg{,s)SI{-s) = I for all s 

2) The M X M matrix 

Sgil - SgSl)-^{I - Sl){I - SgSl)-^S 2 i (13) 

has equal singular values for all s = joo 

3) / — Sl{so)Sg{so) is non-singular 

4) S'J(— s) — S'g(s) has no zeros in the RHP 

where sq is defined in (6). These four conditions correspond to 
four possible impediments to achieving the bounds. Meeting 
all the conditions is sufficient to attain equality in the bounds, 
but Conditions 1, 2 and 4 are also necessary. The N x N 
matrix Sg{s) (see Figure 1) plays a prominent role in these 
conditions and can be readily measured or modeled by con¬ 
necting pairs of output ports of the matching network to a 
network analyzer while terminating its remaining ports with 
characteristic impedances. 

These conditions have physical interpretations. Condition 1 
is satisfied for lossless networks since we are guaranteed that 
<S' 2 i(s)«S'Ji(—s) + 822 ( 3 )Sj 2 {—s) = I for all s because S{s) 
is a para-unitary matrix, and Sg(s) = 822 ( 3 )- 

Condition 2 requires the singular values of (13) to be equal 
for all s = joo- When this is satisfied, the total dissipated 
power at the loads depends only on the total incident power 
|jai(jw)|p, and not the “direction” of di(joo)- 

Condition 3 is a “non-degenerate” condition that we illus¬ 
trate with an example: Suppose the loads are capacitive to 
ground, and hence are reflective with 8 l(oo) = —I- If the 
output impedance of the matching network is also capacitive, 
then 8 g(oo) = —I and Condition 3 is violated. Hence, in 
this example, a matching network that wants to achieve high 
bandwidth should avoid capacitive output impedance. It turns 
out that Condition 3 is superfluous in Bound 2 because 8 l(s) 
and S'g(s) are bounded matrices; hence 8 l(s) 8 g( 3 ) is also 
bounded and / — 8l(so)8g(3o) is always non-singular for 
Re{so} > 0 [24, 7.22]. Additional details on Condition 3 can 
be found in [9]. 

Condition 4 is a “minimum-phase” condition since the RHP 
zeros of 8'£(—s) — S'g(s) are the same as the RHP zeros of 
S'gm(s), which involves a “Darlington equivalent” network 
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representation of the loads [28], [29] as described in Part II. 
However, knowledge of the Darlington equivalent network is 
not needed to check this condition. We have more to say 
about network architectures that cannot meet this condition 
in Section IV-C. 

Assuming Conditions 1^ are met by the matching network, 
we obtain equality in (11), (12) if the network also achieves 
the ideal response r{uj) = r for a; € [oji,oj 2 ] and r(w) = 1 
elsewhere. Any network such that r(w) < 1 for oj ^ [a;i,a; 2 ] 
sustains a non-negative “shaping loss” which is the difference 
between the integral over all uj of the left-hand sides of (7)-(9) 
versus the integral over [wi, a; 2 ]: 

shaping loss = 

r*00 1 pCJ2 1 

/ /(w) log- / /(w)log-^dw. (14) 

Jo r{uj) r{uj) 

For a network whose shaping loss is positive, some perfor¬ 
mance is lost in the band of interest because either a ;2 — Wi 
could potentially be made larger, or t could be made s- 
maller, by redesigning the network to have larger r{uj) for 
oj ^ [a;i,a; 2 ]. We provide an example of the shaping loss 
computation in Part II, Section III-A. 

D. How to obtain Sl{s) and its poles and zeros 

If analytical circuit models for the loads are known, Sl{s) 
is uniquely determined by the standard Laplace transform rep¬ 
resentations of the model impedance or admittance matrix and 
using the formula S'l(s) = (^l(s) + ZqI)~^{Zl{s) — ZqI). 
Examples of this are presented in Section IV-B. 

Absent an analytical model, numerical methods that model 
the Sl{s) of the loads are needed to extract its poles and 
zeros. The modeling methods are discussed and illustrated by 
example in Part II. 

IV. Bandwidth Analysis 

The bounds yield a variety of conclusions when they are 
applied to various system configurations. We first examine 
decoupled loads in Section IV-A and then coupled loads in 
IV-B. We then identify network architectures that can achieve 
the bounds in Section IV-C. 

A. Decoupled loads 

Let Sl{s) = Si{s)I where Si{s) is a scalar and / is an Wx 
N identity matrix. Thus, there are N decoupled identical loads. 
Let Si{s) satisfy (6) for some sq, and let Bi be computed using 
Table I applied to Si{s) for M = 1. Then Sl{s) satisfies (6) 
at the same has poles and zeros at the same locations 

as Si{s), each with multiplicity N. It follows from Table I 
that the bound for N loads is N times the bound for a single 
load. We have proven the following theorem. 

Theorem 1: For N identical decoupled loads driven by M 
sources 

/■“ IN 

/ /(w)log-^dw < —B;. (15) 

Jo rluj) M 

If N = M, (15) is simply Bi. This is not surprising because 
we have assumed the sources and loads are decoupled and 


therefore we are essentially examining N identical isolated 
systems, each with bound Bi. However, (15) scales linearly 
with N for a fixed M. This formula suggests that N identical 
decoupled loads can achieve N times the bandwidth of one 
load as long as the matching network is designed properly. 

Matching network architectures can be thought of in terms 
of their ability to achieve linear-in-W bandwidth behavior 
when used with a set of loads with linear-in-W behavior. 
In Section IV-C we show that certain network architectures 
provably have linear-in-A^ behavior, while others do not. As 
shown in the next section, this linear-in-A^ behavior extends 
to coupled loads under some conditions. 


B. Circulantly-coupled loads 

Loads with circulant Sl{s) are especially easy to manipu¬ 
late mathematically because, while the eigenvalues of Sl{s) 
depend on s, its eigenvectors do not. This makes computing 
PL,i and ZL,i straightforward. Load structures that lead to sym¬ 
metric circulant (s) are identical, and the coupling between 
load i and neighbor j depends on min{|z — j\,N — \i — j\}- 

Two identical loads automatically have circulant symmetry, 
as long as they have reciprocal coupling. Three or more loads 
can be placed in a circular or spherical arrangement to yield 
circulant S'l(s). 

1) Loads exemplifying Bound 1: Figure 2(a) illustrates N 
circulantly-coupled loads. The loads consist of an W-port 
inductive network with impedance matrix sZi and an W-port 
series LC network with impedance matrix ^ — Zic, resonat¬ 

ing at jwQ. These networks are connected in parallel to each 
other, and terminated by isolated characteristic impedances Zq. 
The parallel A^-port networks are shown in Figure 2(b,c); the 
inductive network has each port grounded through Lq and 
every pair of ports i and j is connected through L^, where 
£ = min{|z — j\, N — \i — j\} is the “distance” between ports i 
and j. The series LC network has each port grounded through 
Lq and Cq, and every pair of ports i and j is connected through 
L'g and C'^, where LqCq = y/lyC^ = l/w§. The N x N 
matrices Zi and Zic are circulant symmetric. 

Since circulant matrices have eigenvectors that are columns 
of a discrete Fourier transform (DFT) matrix, we obtain the 
following eigenvalue decompositions: 

Z, = WAiW^, Zi, = WAi,W^, 


where 


W=^ 

Vn 


1 e « 


j2,r(W-l) 

1 e « 


1 

j2,r(W-l) 

e « 


( 16 ) 


_ j27r(W-l)^ 

e « 


is the N X N unitary DFT matrix, and A/ and Aic are real 
positive diagonal matrices representing the eigenvalues of Z/ 
and Zic, respectively. 

The impedance matrix of the loads is 


Zl{s) 



2s 




-zfc" + 
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Fig. 2. (a) N circulantly-coupled loads consisting of an inductive network 

and an LC network in parallel, terminated in a set of characteristic impedances 
Zq. The inductive and the LC portions of the network are shown in (b,c), and 
have impedance matrices sZi and ^ respectively. The loads are 

coupled to ground through Lq and series Lq and Cq, and every load pair is 
coupled through and series and C'^, where i is the “distance” between 
the pair. 


Then Sl{s) can be obtained from Sl{s) = {Zl{s) + 
ZoI)-^{Zl{s) - ZqI), which is 


Sl{s) = W 


_~-^o[(2A^^^a;Q+A^ ^uJq] _ 

2s3+^o(2A^^^cdQ+A^ ^)s2+2ci;QS+ZoAj ^oJq 


w 


H 


(17) 


Note, we write (•)“^ as when we take inverses of diag¬ 
onal matrices. Because Sl{s) is a circulant matrix, only its 
eigenvalues depend on s and the poles and zeros of Sl{s) are 
therefore the poles and zeros of the individual eigenvalues. 

All N loads present short circuits to ground at sg = 0 
and jujo, where the loads become perfect reflectors. Hence, 
Sl{0) = <S'l(7Wo) = —I, and (6) is satisfied. We can then 
obtain two distinct bounds by applying Bound 1 for sg = 0 
and jujQ. Using pL,i, ZL,i calculated from (17), Bound 1 yields 


^OO -| 

.-2 7 ^ 


uj log ■ 


■(u;) 


duJ<^Y' 
- 2M ^ 


— 


TT • trZ; 


MZg 




< 


2M 


E 


2 

— 2A;e,i 
ZgWg MZqOJq 


'(w) 


27r • trZ; 


Ic 


2 ■ 


, (18a) 


(18b) 


The zth diagonal element of Zi represents the inductance 
relative to ground of port i of the inductive network in 
Figure 2(b), measured with the remaining ports open. Because 
Zi is circulant its diagonal elements are all equal. Hence 
(l/iV)trZ/ = Ueq.Af, where Leq,Ar is the inductance of any 
port relative to ground. A similar conclusion holds for the 
LC network in Figure 2(c). Let the equivalent series LC 
branch of any port relative to ground have inductance ^ 
and capacitance C'^^, where C'q.ArC'eq.tv = l/'^o- Then 
(l/V)trZ/c = w^Ueq.AT + Wc Can rewrite (18) as 


1 


oj log . . du) < 
r{uj) 


MZn '■ 




< 


Ntt 


MZi 


2 


r(a;) 


eq.AT I 2r" 

0 \ ^O^eq.AT 


(19a) 


(19b) 


Let Leq,Ar, ^ and C'^ jy approach respective limits Leq, L'^ 
and C'q as N ^ oo. We have proven the following theorem. 

Theorem 2: The linear-in-V behavior shown in Theorem 
1 for decoupled loads extends to circulantly-coupled loads 
provided Leq > 0 and > 0 (or C'^ < oo). 

The conditions Leq > 0 and > 0 (or C'q < oo) are 
equivalent to ensuring that the inductance (or capacitance) of 
any port relative to ground does not go to zero (or infinity) as 
W —> 00, and hence there are not “too many” parallel paths 
to ground from any port. This is equivalent to imposing a 
condition that the coupling between loads not be “too strong”. 

For a given N and M, the bounds in (19) for Figure 2 
increase as inductive and capacitive coupling components are 
removed, since Leq.Af increases and approaches Lq as all cross- 
inductive elements are removed; similarly ^ increases to 
Lq and jy decreases to Cq because Leq.w, ^eq,N 

are obtained as Lq, Lq and Cq in parallel with the remaining 
parts of the networks. 

However, we cannot conclude that coupling always has a 
negative effect on bandwidth. In fact. Part II of this paper 
shows that coupling between dipole antennas has a non¬ 
monotonic effect on the bandwidth bound. This dichotomy 
in behavior between the model in Figure 2 and the dipoles in 
Part II is not contradictory because there is no requirement 
that dipoles should be modeled by fixed Lq, Lq, and Cq as 
the coupling between them changes. 

2) Loads exemplifying Bound 2: Figure 3 illustrates N 
resistors Zq terminated with two parallel A^-port networks: 
one is resistive with N x N circulant impedance matrix Z^, 
and the other is capacitive with N x N circulant impedance 
matrix ^Z^- The impedance matrix of the loads Zi,{s) is 

Zl{s) = ZqI + {Zj. ^ + sZ^ 


Let Ar and be the eigenvalue matrices of Zj. and Z^,. Then 
Sl{s) is 

Ar-A,- 


Sl{s) ^2ZQjy^s -f 2ZoAc -f KK 

where W is given in (16). 

We note that 

I 


■W 


H 


( 20 ) 


. -4Zq A^S^ -I- 4Zg(Zg/ -f Ar)A^ „ 


AcCZo^ZoI+Af) j 

Z„A.. - = O-of 


lU". 


Let the component values satisfy 


for 


some tTg > 0. Then (6) is satisfied for sg = cg. We substitute 

Ac = ^zllzlT+K^) 

^ ° log -^dio 


(Jq +0J 


riuj) 
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Fig. 3. (a) circulantly-coupled loads consisting of N resistors Zq in series 
with a pai'allel resistive network Zr and a parallel capacitive network ^Zc- 
The resistive and the capacitive portions of the network are shown in (b,c), 
and have respective impedance matrices Zr and The loads are coupled 
to ground through Rq and Co, and every load pair is coupled through 
and C^, where i is the “distance” between the pair. 





Fig. 4. An RF system where one source drives N loads having S- 
matrix Sl through a constrained matching network consisting of a prescribed 
(Af+l)-port network Sq an arbitrary network S\ (the complex-frequency 
argument s is omitted). We use to denote the S-parameter seen from 

the input of So- 


the moduli of the eigenvalues of a real symmetric matrix are 
equal to its singular values, at least one eigenvalue satishes 
|AG(ja;)| = 1. We conclude that reciprocal broadband splitters 
or couplers that yield real symmetric Sc{s) cannot be used 
to achieve equality in (15) when M < N and the loads are 
decoupled. We see an example of this in the next section. 


— TT 

< ivTT iog det — 
“ 2M ^ 2Zo 


TT 

W 


N 


^log 


2=1 


Ar + 2-y/ ^o(-^0 + Ar) 
2Zq -f + 2^Zq(Zq + Ar,i) 

■ 

Z\r.r 


( 21 ) 


where A^ ^ are the diagonal elements of A^-. The linear-in-V 
behavior of these loads depends on A^.i,..., A^^at; we do not 
pursue this analysis any further here. 


C. Ability of network architecture to achieve bounds 

Matching network architectures that cannot achieve equality 
in the bounds should be identihed whenever possible and 
removed from consideration in large-bandwidth applications. 
The following theorem shows how Sg{s) dehned in Figure 1 
can be compared with Sl{s) to identify such networks. We 
thereby make Condition 4 in Section III-C more physically 
tangible. 

Theorem 3: For matching networks satisfying Condition 3, 
if Sg{s) has an eigenvector in common with Sl{s), and the 
associated eigenvalue Ag(s) satishes |AG(ju;)| = 1 for all ui, 
then Condition 4 cannot be satished and therefore the bounds 
cannot be achieved. 

Proof: See the Appendix. ■ 

Theorem 3 can be applied to the example of decoupled loads 
in Section IV-A. We assume \ < M < N and Condition 3 is 
satished. Since Sl{s) = Si{s)I, any vector is an eigenvector 
of S'l(s). Theorem 3 then implies that any S'g(s) that has 
|Ag(jw)| = 1 cannot achieve equality in (15). 

For example, networks where Scis) is a real symmetric 
matrix cannot achieve equality. This follows because even if 
the network is lossless and Condition 1 is satished, implying 
that S' 2 i(s)S'Ji(—s) + «S'g(s)5'q(— s) = I for all s, then 
because S' 2 i(s) is an V x M matrix and Sg{s) is an N x N 
matrix, Sg{s) has at least N — M unit singular values. Since 


V. One Source AND Multiple Loads 

We probe the special case of one source more deeply 
with some examples in this section. We analyze constrained 
networks in Section V-A. Decoupled and coupled loads are 
examined in Sections V-B and V-C, leading to a class of non¬ 
reciprocal networks called “determinant” networks. 

A. Constrained matching network 

Let a portion of the network be constrained or prescribed 
to have a particular structure. For example, the constrained 
portion might include ideal circulators or power splitters. 
The prescribed portion and loads then together constitute an 
“equivalent load” and the achievable bandwidth is determined 
by the characteristics of this load. 

Figure 4 illustrates an example of such a structure for M = 
1, where a prescribed (V+ l)-port network is combined with 
the loads and forms an equivalent one-port load. Let the S- 
matrix of the prescribed network be 

1 N 

So(s) = ^ ( ‘^ 041 ( 5 ) -S'o, 12 (s) j 

w V 5-0,21 (s) 5 o,22(s);' 

Then the S-parameter of the equivalent load is 
5L,eq(s) = 5o,ll(s) 

+ 5o.i2(s)5l(s)(/- 5o.22(s)5l(s))-15o.2i(s). (22) 

The remaining unspecihed portion of the network 5i(s) is 
connected to the input port of SL,eq{s), and we may ask what 
bandwidth is attainable by 5i(s). 

Generally, we cannot apply the bounds directly to SL,eq{s)- 
Unlike 5l(s) where power is either dissipated or reflected 
by the loads, power in 5L,eq(s) can also be absorbed by the 
constrained portion of the network when 5o(s) is lossy. Let 
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ri{uj) denote the ratio between the power dissipated by S'l(s) 
and the power delivered to S'L,eq(s), defined as 

\\a2iju;)r-\Mju;)r 

^ McW-IMjcoW^ 

where ( 02 (jw), b 2 (juj)) and (aoijuj), bo(juj)) are the (inci¬ 
dent, reflected) signals from Sl(s) and 5'L,eq(s), respectively. 
Then 0 < r/(uj) < 1 and 

r](uj) = 

- So, 22 SLr"(I - S^Sl)(I - So, 22 Sl)-^So, 21 

1-^^eq^L.eq 

(23) 


which depends only on the prescribed network and the loads. 

We assume S'L,eq(s) satisfies (6) for some sq; this does 
not require Sl(s) to satisfy (6) for the same sq- Then (10) 
becomes the following bound. 

Theorem 4: For loads matched by a constrained network, 

/ /(W) log j - -duj < Beq (24) 

Jo r^{uj)+r){uj)-1 

where Bgq is the right-hand side of (10) applied to S'L,eq(s)- 
Proof: Using Figure 4, we define 

2, . , 11000*^)11" - l| 4 o(i‘^)ll" 

= ‘-WiO#— 

as the ratio between power not delivered to S'L_eq(s) and the 
incident power from the source. Then it follows that 

2 . 2 , \\d2{3u:)r-\\b2{j^)r 

_ \\ao{ju:)r-\\boiju;)r ||« 2 (j^)f - 

||ai(jo;)P WaoijcoW-WboijivW 

= 1 - (1 - 


Equivalently, we have req(a;) = 

We apply (10) to the equivalent load SL,gq{s). The result is 
an inequality on req(w): 

[ /(w) log - ]-^duj < Beq, 

Jo ^eq(w) 

where i3eq is the right-hand side of (10) applied to SL,eqis). 
Replacing req(a;) with \J gives us (24). ■ 

We use beamforming as an example application. Let the 
desired amplitude and phase relationship of the antennas be 
denoted by the iV x 1 unit real-rational vector v{s) for s = joj. 
Then 


So{s) 


0 bF{s) 

v{s) 0 


(25) 


denotes the S-matrix of a reciprocal one-to-iV power divider 
that constrains the incident signal 02 (jw) to be aligned 
with v{juj). We are not concerned with the so-called gain- 
bandwidth product of a phased array [30] which measures its 
ability to maintain far-field gain across a range of frequencies. 


Theorem 5: Let the beamforming vector v{s) be a real 
constant unit eigenvector of Sl{s). Then 

/ (W) log -/-i- duj < Beq, (26) 

r(a;) 

where Bgq is calculated using S'L.eq(s) = u^(s)5'l(s)F(s). 

Proof: We substitute (25) into (22) and (23), where 
'S' 0,11 = 0, 5'o,12 = (s), <S'o ,21 = v{s) <S'o ,22 = 0. 

Then (22) yields 

-S'L.eq(s) = v'^(s)S'l(s)F(s), 



and (23) yields 

- S^Sifv _ v^v- v^S^Slv 
1 - 5'l^eq'S'i.eq 1 - Sj^V* iFS lv' 


Because v{s) is a real unit eigenvector of Sl{s), it is readily 
verified that ri{uj) = 1. The result then follows by applying 
Theorem 4. ■ 

Section III-B in Part II of this paper provides an example 
where Bgq in (26) varies with the choice of v{s). 

We note that r]{uj) = 1 is obtained for any lossless So{s), 
in which case (26) also applies. When r]{uj) = 1 and Sl{s) 
satisfies (6) for some sq, then 5'i,_eq(s) also satisfies (6) for 
the same sq; the converse is not true. 

Compared with (10), (26) is smaller for N > 1. To see 
this more explicitly, let the loads be decoupled as in Section 
IV-A, whence S'l(s) = Si{s)I. Then any real constant unit 
vector v{s) satisfies r]{uj) = 1, and the equivalent load satisfies 
Sb.eqis) = Si{s). Then Theorem 5 yields 


f{<^)^og-^duj < Bi, 
r(a;) 


(27) 


where Bi is the right-hand side of (10) applied to Si{s) for 
M = 1, Clearly, Bi < NBi, which is the bound obtained in 
(15) for M = 1. We conclude that linear-in-A^ behavior is not 
generally achieved for the beamforming structure (25). This 
conclusion is consistent with Theorem 3. 


B. Decoupled loads 

Section IV-C says that full bandwidth cannot be achieved 
when reciprocal broadband splitters are used to drive > 1 
loads when M = 1. We instead consider a matching network 
consisting of non-reciprocal couplers in Ligure 5(a) where an 
{N + l)-port network comprising N — 1 circulators as the 
constrained part of the network is displayed. The resulting 
{N + 1) X {N + 1) S-matrix is 


^o(s) 


0 0 ••• 0 1 

1 0 ••• 0 0 

01 0 ••• 0 

: 0 : 
0 ••• 0 10 


The lower right N x N block of this matrix, which cor¬ 
responds to Sa{s), has N zero eigenvalues and therefore 
does not satisfy the conditions of Theorem 3. Let S'L,eq(s) 
be the S-matrix of the equivalent one-port load seen at the 
input of the circulators. Because the circulators are lossless. 
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Sl{s) into det5'i(s) a “determinant network”. Determinant 
networks are linear-in-A^. 

Let W be defined as in (16) and have columns wi,... W]s[. 
Define Wi = [w 2 W 3 ■ ■ ■ wn and the {N + 1) x {N + 1) 
matrix 


So{s) = 


0 

Wi 


W§ 


(28) 


Fig. 5. (a) — 1 circulators are used to match N decoupled loads to 

M = 1 source, and achieve linear-in-behavior. We use 5'L,eq(-5) to denote 
the S-parameter of the equivalent one-port load seen at the input of the 
circulators, (b) N decoupled loads ai‘e driven by M = 1 source through 
N two-port networks and TV — 1 broadband circulators. Each of the two-ports 
matches 1/TV of the total bandwidth; network i passes band i and reflects the 
remaining portions. The circulators combine all TV passbands and achieve a 
total bandwidth TV times of the bandwidth of each two-port network. 


r]{uj) = 1 in this constrained network. It is readily calculated 
that S'L,eq(s) = which has poles and zeros at the 

same locations as Si{s), each with multiplicity N. Therefore, 
according to (10), the bandwidth achievable with this equiv¬ 
alent S'L.eq(s) is N times that achievable by Si{s). Hence, 
linear-in-V network behavior for N decoupled loads can be 
achieved by ideal circulators. 

This circulator-based network architecture presents “multi¬ 
ple opportunities” for the energy that is refiected from any 
one load to be forwarded to the next for another attempt at 
transmission. We are ignoring the insertion losses associated 
with cascading circulators in such an arrangement. 

Figure 5(b) has a structure similar to Figure 5(a) but matches 
1/iV of the total bandwidth in an orderly fashion. The first 
network passes the lowest portion of the band and refiects 
all the remaining portions, which are passed to the next 
circulator which matches the next portion, and so on. This 
structure resembles a channelizer [31] in the sense that there 
are multiple output ports where each port contains a portion of 
the total bandwidth. Although we have drawn these networks 
as having near-ideal fiat frequency responses, this aspect is not 
crucial. An explicit example of this type of network for a pair 
of dipoles is presented in Section III-A in Part II. 


Notice that Wi is missing the column wi. Then Sq{s) is 
constant and lossless. The following theorem says that it is 
also a determinant network. 

Theorem 6 : Let the network described by (28) have its N 
outputs connected to any circulantly-coupled set of loads with 
S-matrix Sl{s). Then its input has equivalent S-parameter 
5'L,eq(s) = det S'l(s). 

Proof: We let L = W^Wi and use the eigenval¬ 
ue decomposition Sl{s) = WA.{s)W^ where A(s) = 
diag(Ai(s),..., Ajv(s)) is a diagonal matrix of eigenvalues. 
Then (22) yields 

SlM = w%Sl{,s){I - WiW^S l{s))-^wi 

= w^WK{s)W^{I - WLK{s)W^)-^wi 
= [0 ••• 0 \n{s)]{I-LK{s))-^ 


where 


I-LK{s) = 


<[10 

1 

-Ai(s) 

0 

0 


0 


0 

1 

-A2(s) 


0 —Ajv-i(s) 1 


It follows that SL,eq{s) is AAr(s) times the (W, 1) entry of 
(/ - LA(s))“^ This gives SL,eqis) = HTILi An(s). ■ 

An intuitive explanation of the operation of the network 
is as follows. The network (28) orients the power from the 
source along the first eigenvector wi. Energy is then refiected 
by the loads with amplitude Ai(s) along wq, at which point 
the determinant network refiects it entirely back to the loads, 
but with orientation W 2 - This is refiected by the loads with 
amplitude A 2 (s) which is returned by the network to the 
loads reoriented along ws, and so on. The last eigenvector 
wn refiected by the loads is returned to the source. The result 
is an overall S-parameter with value 11^=1 '''n(s), which is 
the determinant of S'l(s). This description also explains why 
a determinant network is not unique. 

For example, for N = 2, 


C. Circulantly-coupled loads: Determinant networks 

The previous section gave a circulator-based architecture 
for achieving linear-in-V bandwidth for N uncoupled loads. 
We now show how to handle circulantly coupled loads such 
as examined in Section IV-B, by demonstrating a network 
architecture that converts the multiport system with S-matrix 
Sl{s) into a single-port system with S-parameter SL,eq{s) = 
detS' l(s). Since detiS'L(s) has the same poles and zeros as 
S'l(s) if S'l(s) has no cancelling poles and zeros, SL,eq{s) has 
the same bound as Sl{s). We denote any network that converts 


So{s) 


0 

1 

f 


1 

2 


1 

2 



(29) 


As can by seen by this S-matrix, the network works by first 
exciting the even mode of the coupled loads. Any refiected 
power, which is also in an even mode, is then converted into 
an odd-mode excitation. 

One possible realization of (29) is given in Figure 6, which 
utilizes an ideal transformer and gyrator. A gyrator with 
resistance Zq is transparent to forward propagating waves but 













































11 



Fig. 6. Realization of the network (29) that includes a transformer and 
non-reciprocal gyrator connected to two coupled antennas. 


behaves as a tt phase shift for reverse waves. The network 
in Figure 6 converts the wave incident on the transformer 
primary into an even-mode excitation for the loads on its 
secondaries. Any power reflected by the loads, which is also in 
an even mode, is converted into an odd mode by the gyrator. 
The secondaries of the transformer acts like an open circuit 
to this odd-mode signal. The signal is therefore reflected by 
the transformer back to the loads as an odd-mode excitation. 
Finally, any odd-mode excitation reflected by the antennas is 
then converted by the gyrator into an even-mode signal and 
passed by the transformer back to the source. 

The network in Figure 6 can achieve the bound in (10) for 
a single source and two identical loads that are reciprocally 
coupled. This network would also work well for the decoupled 
loads in Section IV-A since decoupled loads are trivially recip¬ 
rocal. Note that the circulator structure proposed in Figure 5 
is a determinant network when the loads are decoupled but is 
otherwise not. Hence, the network in Figure 6 generalizes the 
networks in Figure 5 for N = 2. 

VI. Conclusions 

We have presented bandwidth analyses for multiport RF 
systems using bandwidth upper bounds that apply to an 
arbitrary number of sources and loads, and allows arbitrary 
coupling between the loads. Conditions for achieving equality 
in the bounds were discussed. We demonstrated that the 
bandwidth bounds scale generally as N/M for M sources 
and N loads. We focused on loads whose scattering matrices 
can be expressed analytically with rational functions. The case 
of one source and many loads was examined deeply. 

In Part 11 of this paper, we apply the upper bounds to 
realistic loads whose scattering parameters are expressed nu¬ 
merically and therefore need to be approximated by rational 
functions. The accuracy needed of such approximations is 
analyzed. Some of the effects of coupling are examined 
in detail, and complete derivations of the bounds are also 
provided. 

Appendix 

Proof of Theorem 3 

Let v{s) be the eigenvector of Sl{s) and Sc{s) associated 
with eigenvalues Xl{s) and Ag(s). Because Sl{s), Sc{s) 
are bounded matrices, Al(s),Ag(s) are bounded functions. 
Since |Ag(ju;)| = 1, Ag(s) is an all-pass function with all 
poles in the LHP and zeros in the RHP. Moreover, because 
I — S'£{—s)Sl{s) is full normal rank and Sl{s) satisfies 


(6), Al(s) is non-constant and satisfies Al(— so)Ai(so) = 1 
where sg is defined in (6). Thus, |Al(s)| < 1 for Re{s} > 0, 
|Ag(s)| > 1 for Re{s} < 0 and |Ag(s)| < 1 for Re{s} > 0; 
the equalities holds if Ag(s) is a constant. 

We start by showing {juj)SG{joj) = XG{joj)v^{joj). 
Since Xcijoj) and v{joj) are an eigenvalue-eigenvector pair 
of Scijuj), 

'y^(jw)(5'G(jw) - XGijuj)I)v{juj) = 0. 

Suppose u^iju}) = v^{juj){SGiji^) - XgUuj)!) ^ 0, then 
u{juj)J-v{juj), and 

\\v^{ju})SG{juj)\\ = ||AG(jw)F^(ja;) -f u^(ja;)|| 

> l|AG(jw)F^(ja;)|| = ||i;^(ja;)||. 

But because Sg{s) is bounded, all singular values of S'g(jw) 
are no greater than one. This contradicts the inequality 
above. Hence, u^{juj) = 0, and we have {jijj)SG{joj) = 
XG{juj)v^{juj). 

Let F'(s) be the extension of v*{juj) to the WCP We then 
show that Condition 4 cannot be achieved when Condition 
3 is satisfied. We need to show {s){Sj^{—s) — SGis)) = 
(Al(— s) —AG(s))f^^(s) = 0 for some s in the RHP. It suffices 
to show Xl{—s) — Xg{s) = 0. 

Since |Al(— s)| < 1 and |Ag(s)| > 1 for Re{s} < 0, all 
zeros of Al(— s) — Ag(s) = 0 must locate in Re{s} > 0. 
Both Al(s) and Ag(s) are rational functions, so we let their 
degrees be ni and n 2 , respectively. Since all the poles of 
Xl{s) and Ag(s) are in the LHP, the poles of Xl{—s) and 
Ag(s) do not coincide. Hence Xl{—s) — Ag(s) has degree 
Til + Moreover, |Ai(jw)| = 1 has at most ni solutions; 
this is because 1 — Al(—s)Al(s) has degree at most 2ni, 
and any zeros of 1 — Al(—s)Al(s) on the imaginary axis 
have multiplicity at least two. Since |AG(ju;)| = 1 for all oj, 
Xl{—s) — Ag(s) has at most ni zeros on the imaginary axis. 
This leaves at least n 2 zeros in the RHP, and we have proven 
the theorem when n 2 > 1. 

When n 2 = 0, we prove by contradiction. Suppose all ni 
zeros of Xl{—s) — Ag{s) are on the imaginary axis. Then the 
zeros of 1— Al(—s)Al(s) are at the same locations as the zeros 
of Xl{—s) — Xg{s). Because we assume Al(—so)Al(so) = 1 
for some Sg must be one of the zeros of Xl{—s) — Ag(s). 
Thus, we have Ag(so) = Xl{-so), and 1 - Ag(so)Al(so) = 
0. This contradicts Condition 3. Hence, Al(— s) — Ag(s) has 
zeros in the RHP when n 2 = 0. Thus, Condition 4 is not met, 
and since this condition is necessary for achieving the bounds, 
the bounds cannot be met with equality. ■ 
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Bandwidth Analysis of Multiport Radio-Frequency 

Systems—Part II 

Ding Nie, Member, IEEE, and Bertrand M. Hochwald, Eellow, IEEE 


Abstract —We analyze the bandwidth over which power can 
be transferred from multiple radio-frequency sources to multiple 
loads through a passive multiport matching network. This is the 
second part of a two-part paper. In the first part we introduce 
broadband multiport matching upper bounds and apply them to 
determine the bandwidth of loads and network structures whose 
scattering parameters can be expressed analytically with rational 
functions. In this second part, we apply the bounds to loads, 
such as antennas, whose scattering parameters are obtained by 
measurement or simulation. We focus on the effects of coupling 
on bandwidth. Since the bounds require frequency responses 
that are rational functions, we provide guidelines on how to 
obtain rational approximations for arbitrary loads. Complete 
derivations of the bounds are also provided. 

Index Terms —Bandwidth, Bode-Fano bounds, broadband 
matching bounds, non-reciprocal networks, passive matching 
networks, radio-frequency coupling 


I. Introduction 

We analyze the bandwidth over which power can be trans¬ 
ferred from M sources through a passive matching network 
to N loads in radio-frequency (RF) systems. Our principal 
tools include broadband matching upper bounds, which are 
presented in Part I of this paper [1], and briefly repeated herein 
for easy reference. These bounds extend the classical Bode- 
Fano results [2], [3], which apply to a single source and load. 
The bounds depend only on the scattering matrix (S-matrix) 
of the loads. The loads may be electromagnetically coupled to 
each other. 

In this part, we apply the bounds to loads whose S-matrix is 
expressed numerically, found through either measurements or 
simulations. The bounds require rational models of the loads 
that can be obtained through fitting of the numerical data. 
We provide guidelines for assessing the accuracy of a rational 
model. 

Of particular interest is the effect of coupling on the 
bandwidth of multiple antennas. One example demonstrates 
that the bandwidth bound of a pair of dipoles fluctuates non- 
monotonically with the distance between them, and, with the 
right amount of coupling, can be significantly greater than for 
decoupled dipoles. We focus on the “input bandwidth” of the 
loads, and are not concerned with how efficiently the loads 
use their input power. 
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kirknie@gmail.com). 

B. M. Hochwald is with the Department of Electrical Engineering, 
University of Notre Dame, Notre Dame, IN 46556, USA (e-mail: b- 
hochwald@nd.edu). 
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We also include complete derivations of the bounds. These 
derivations are not needed to apply the bounds but are of use 
in understanding how the conditions for equality in the bound 
apply. 

The notations and definitions are identical to those in 
Part I. We assume M decoupled sources with characteristic 
impedance Zq drive N loads that have a frequency-dependent 
N X N S-matrix S'l(s), where s = a + juj is the extension 
of frequency joj to the whole complex plane (WCP). An 
(M + 7V)-port passive matching network is inserted between 
the sources and loads. Our measure of the quality of matching 
is the power loss ratio G [0,1], defined as the ratio 

between the total power lost, including insertion and return 
losses, and the total incident power from the sources at 
frequency uj. The insertion loss is the power dissipated in the 
matching network; the return loss is the power reflected to the 
sources. Since the loads are potentially coupled, the quantity 
r^(w) captures the possibility that incident power from one 
source can reflect back to itself or to other sources. We then 
define the bandwidth wbw(t, as the frequency range for 
which r{uj) (the positive square root of is no greater 

than a threshold r > 0 in the vicinity of a design frequency 
ujd. The mathematical definitions of r(a;) and wbw(t, can 
be found in Part I. 

We generally employ the properties of S-matrices of passive 
real networks to prove the bounds; these properties include 
real-rational, Hurwitzian and bounded. We also use the defi¬ 
nitions of poles and zeros of rational matrices. We use LHP 
to denote the left-half complex plane (Re{s} < 0) and RHP 
to denote the right-half (Re{s} >0). 

A. Summary of bounds 

There are three bounds which are summarized as follows. 
If 

SI{-so)Sl{so) = I (1) 

is satisfied for some Re{so} > 0, then the following inequality 
holds for any passive matching network: 

/(w)log-^da; < (2) 

r(a;) 

where /(w) and B depend on the location of sq and are given 
in Table I. To compute the bound B, the poles and zeros of 
Sl{s), denoted as PL,i, ZL,i, are also needed. We refer to rows 
1, 2, and 3 in Table I as Bounds 1, 2, and 3. 

The bounds depend on the poles and zeros of the real- 
rational matrix Sl{s), which we explain how to obtain in 
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TABLE I 

Forms of f{juj) and B in (2) for different locations of so, 

WHERE ARE THE POLES AND ZEROS OF Sl{s). 


so = io20 

/(co) = |[(a 2 o - co) 2 + 2] 

B [Ei(PL,i -io2o)“l + Ei(2L,i + 2 ^ 00 )“!] 

Re{so} > 0 

/(a;) = iRe[(so - juj) ^ + (sq + juj) 

wlog|det5L(so)-[{;|::!:^;:j 


SO = 00 

/( 02 ) = 1 

i PL,i + Xli ^L,i) 



Section II. We focus on the steps needed to fit realistic loads 
with rational models. Bandwidth analyses of various loads are 
then presented in Section III. The bounds are corollaries of 
Theorems 2 and 3, which are presented in Section IV. We 
conclude in Section V. 

II. Obtaining Sl{s) and its Poles and Zeros 

The matrix S'l(s) = {Zl{s) + ZoI)~^{Zl{s) - ZqI) is 

uniquely determined if the impedance matrix of the loads 
Zl{s) is known analytically. For many realistic loads, this 
is generally impossible. Instead, either ZL{juj) or SL{juj) 
is often obtained as the result of a numerical simulation or 
measurement, where oj is contained to some frequency range 
of interest. The next section gives an example of how to use 
the results of simulation to find S'l(s). 

A. Finding Sl{s) from numerical simulations of loads 

The following steps can be used to find Sl{s): 

a) Measure or simulate the loads in the frequency range of 
interest [ 021 , 022 ]. Denote the measured response 

b) Find a passive rational S'l(s) such that Sp{joj) is “close 
enough” to S'j^{ju) for u G [ 021 , 022 ]. 

Step a) can be done with standard modeling software such as 
Ansys HFSS in the case of simulations, or a network analyzer 
in the case of measurements. Step b) can be accomplished by 
fitting rational functions to the individual entries of S'^{juj) 
using, for instance, the Matrix Fitting Toolbox [4]-[8] in 
MATLAB. We have more to say about this fitting in the next 
section. 

To compute the right-hand side of (2): 

c) Find an sq where (1) is satisfied for the computed Sl{s). 

d) Calculate the poles and zeros PL,i,ZL,i of Sl{s). 

Step c) requires us to solve S'£{—so)Sl{so) = I. For 
So = ju 2 o, we can instead solve | det S'i,(ju 2 o)| = 1. For Step 
d)^ PL,i,ZL,i can be obtained from the definitions of poles 
and zeros of matrices. In some cases the poles and zeros of 
S'i(s) coincide with the poles and zeros of detS'L(s). One 
way of checking this is presented in [9], which we do not 
repeat here. Then pL,i, ZL,i can be obtained by applying root¬ 
finding algorithms to l/detS'L(s) = 0 and detS'L(s) = 0, 
respectively. The examples shown in Sections III-B-III-E use 
this approach. 

Multiple rational models may exist within the error toler¬ 
ance of for 02 G [ 021 , 022 ]. Each model could satisfy 


(1) for different sq, and generate different bounds. When these 
models yield consistent results, one can build confidence that 
the bounds and models are physically meaningful. When these 
models contradict each other, further investigation is needed 
to determine the source of the inconsistency. An example of 
multiple models is shown in Section III-C. 


B. Analysis of distributed-element loads using rational models 

Since S'l(s) is a real-rational matrix, accurate modeling 
of lumped-circuit loads is straightforward and can be done 
without error. However, many loads such as antennas, trans¬ 
mission lines, and other distributed-element systems include 
time-delays, stubs, and other structures that are often modeled 
using non-rational functions of s. A real-rational S'l(s) is then 
needed that “approximates” the loads with sufficient accuracy 
to obtain a meaningful bandwidth bound. 

The numerical and analytical approximation of distributed- 
element circuits by lumped circuits over a fixed bandwidth is a 
topic that has been studied in many contexts. An early example 
is [10]. A summary of some practical techniques can be found 
in [11]. Recent advances in antenna modeling include [12]. 
Other fields such as control engineering use rational fitting 
extensively [13]. 

The mathematical problem of rational approximation and 
modeling is addressed by Runge’s theorem in complex analy¬ 
sis [14], which states that any analytical function on a subset of 
the complex plane can be fit arbitrarily precisely with rational 
functions. Eor example, the Fade approximation [15], [16] is 
often used to model time delays using rational functions over 
a frequency band of interest. 

We are interested in obtaining an upper bound on 

r^2 1 

/ /(w)log^^dw, (3) 

where rfuj) is the loss experienced when the matching 
network is used with the actual loads S'j^{joj). The process 
outlined in Section II-A creates a rational approximation 
SlU^) for w G [ 021 , 022 ]. Eor an arbitrary matching network, 
(2) gives 

/ /(o 2 )log-^(i 02 < / /(a 2 ) log-^(102 < B, (4) 

r{uj) Jo r{uj) 

where B is computed from S'l(s) and r{uj) is the power loss 
ratio experienced with S'l(s). The difference between the left- 
hand sides of (3) and (4) is 



/(w)log^^do2, 

r'(o 2 ) 


and represents the error in the integral introduced by the 
approximation of Then choosing Sl{s) such that 


r f(uj)\og^du;<6B 
rfuj) 

for some desired tolerance 6B > 0 yields the bound 


(5) 



/(w) log 


———duj < B SB. 
r{uj) 


( 6 ) 



3 


The following theorem provides a first-order approximation 
of 6B for 5 ^( 5 ) that is close to in the frequency band 

of interest. 

Theorem 1: Let SShijoj) = SlUoj) — and define 


p{u) 


1 - 



- 0'L,min(^) 

(1 - 


(7) 


where as^maxi'-^) is the maximum singular value of 5SL{joj), 
and 0 < min('^) < '^L,max(^) < 1 are the minimum and the 

maximum singular values of S'^{jLS), respectively. Then 


5B 





1 - r'^juj) 

r'2(a;) 


p{uj) ) doj 


( 8 ) 


as ||55'z,(jw)||F ^ 0 for w e [wi,W 2 ]. 

Proof: See Appendix A. ■ 

In (8), /(w) is given in Table I, and p(w) is determined 
by S'j^{juj) and 6SL{joj). The notation || • ||f refers to the 
Frobenius norm (sum of squared magnitude of entries) of a 
matrix. In the limit of a perfect model when 6SL{joo) = 0, 
then p{oj) = 0 and SB = 0. In general, SB depends on r'{uj), 
which depends on the matching network. To eliminate the 
dependence of SB on r'{oj), we choose the desired goal of 
r'{uj) « T for w G [a;i,W 2 ], whence 


SB 





p{uj). ) duj. 


(9) 


Note that even if r'{oj) does not achieve t for w G [u;i,W 2 ], 
(9) is still a useful quantity because it upper-bounds SB. 

The steps for evaluating goodness-of-fit are then; 

e) Compute SB in (9) numerically. 

f) Evaluate if the desired tolerance on SB/B is met. 

If the desired tolerance is not met, a smaller SSL{juj) is 
needed. Finally, 

g) Compute the bound (6). 

Theorem 1 can also be used to examine cases where 
the rational fit may fail to generate an accurate bound. For 
example, if <j'i^ ~ 1 in the neighborhood of some w, 

meaning is very reflective, (7) indicates that piuf) 

becomes large. Any attempt to make rfuj) small around this w 
could then potentially make SB in (8) large. Thus, attempting 
to match a load where it is naturally very reflective may lead 
to a loose bound. 


C. Hazard of over-fitting 

In attempting to make SB/B small, care should be taken not 
to “over-fit” the loads. Over-fitting or over-modeling occurs 
when the degrees of the polynomials in the numerator and 
denominator of the rational functions in 5 'f(s) are larger than 
needed to obtain the desired accuracy. As the polynomial 
degrees are increased, the relative error SB/B is decreased 
since SB is decreased. But there is also the hazard that B 
is made unnecessarily large because of the excessive number 
of poles and zeros in S'l(s). As a result, (6) becomes loose 
and Condition 4 for achieving B cannot be satisfied. There is 
a tension between making the rational model accurate while 
still having a small number of poles and zeros. 


One approach to determine if the system is over-fitted is to 
examine the poles and zeros of the entries of Sl{s). Any poles 
and zeros that nearly cancel can have a small effect on SB but 
a large effect on B-, they should be removed and SB/B re¬ 
checked. Simpler models yield tighter bounds, and thus there is 
an incentive to find the minimal model that adequately captures 
the frequency response of the loads in the band of interest. 
A detailed discussion of modeling poles is found in [17]. A 
discussion of the sensitivity of poles and zeros to perturbations 
of the data is found in [18]. 

In the remainder, we use the notation r{uj) when we are 
considering the rational model, and we use r'(w) when we 
are considering the actual loads. In Section III-D, we present 
an example where time delay in a circuit is modeled using 
rational functions. Section III-E computes (6) for four realistic 
coupled antennas. 

III. Bandwidth Analysis of Antennas 

We now calculate broadband matching bounds for five 
examples, all involving amplifiers driving antennas. Aspects of 
coupling and rational approximation are emphasized. The first 
three examples assume that the rational models are accurate 
and therefore compute only (2), while the last two examples 
incorporate the effects of rational fitting and therefore compute 
(6). We begin with two identical decoupled dipoles whose 
rational models are taken from [9]. 


A. Two decoupled dipoles 

We design a broadband network of the type in Figure 5(b) 
in Part I for two decoupled dipoles (N = 2) with one source 
(M = 1). From Theorem 1 in Part I, the bound for two 
decoupled dipoles is twice as large as for a single dipole when 
there is only one source. The geometry of a dipole is shown in 
Figure 1(a), which is half-wavelength at 2.4 GHz. In [9], the 
same dipole is simulated in the range 1-5 GHz, and is modeled 
using an S-parameter model Sfis) normalized by characteristic 
impedance Zg = 50 O that satisfies Sfioo) = —1. We apply 
the model here, and write the S-matrix of two such dipoles as 
Sl{s) = diag(S'i(s),5'i(s)). 

Clearly, Sl{oo) = —I, and (1) is satisfied at sg = oo. The 
poles and zeros of Sfis) are listed in Table III in [9], and 
Bound 3, using (2), for M = 1 gives 

nUJ2 2 \ 

/ log , ' do; < / log , , dcu < 4.93 X 10^°, (10) 

rluj) Jo r{uj) 


where uji = 2 tt x 10® and a ;2 = IOtt x 10® are the lower and 
upper modeling frequencies. 

We are interested in matching the decoupled dipoles to one 
source in the 2-A GHz range. The structure in Figure 2 is 
employed, where Matching Network 1 is tuned for the 2-3 
GHz band and Matching Network 2 is tuned for the 3^ GHz 
band and one circulator is used. These two-port networks are 
designed using the Real Frequency Technique [19], followed 
by a realization using Darlington synthesis [20]. In Figure 1(b), 
the dashed and dash-dot curves show r{uj) of each network 
when connected by itself to a dipole (typically through a balun 
that is not shown). The r(a;) of the entire system, including 
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(a) (b) 

Fig. 1. (a) Geometry of a dipole that is half-wavelength at 2.4 GHz. (b) Two 
decoupled dipoles are matched by the structure shown in Figure 2. Matching 
Network 1 handles 2-3 GHz and Network 2 handles 3^ GHz; their resulting 
r(aj)’s are given in the blue dashed and red dash-dot curves, respectively. The 
overall system r(a;) is shown by the solid curve. 


circulator, is the product of the r(a;) achieved by each, and 
is shown as the solid curve. We seek a power reflection ratio 
of r{uj) < 0.2 (—14 dB). From the figure, we see that the 
bandwidth achievable is u;bw( 0.2, 3 GHz) = 2.35 GHz. The 
achieved integral is 

r^2 1 

/ log--dw = 4.58 X 10^°. (11) 

JuJ^ r{uj) 

There is only a 3.51 x 10® gap between (10) and (11), which 
is due entirely to shaping loss defined in Part I, Section III-C. 
We therefore achieve excellent performance for one source 
and two dipoles with this non-reciprocal network. According 
to Theorem 3 in Part I, since M < N, using a reciprocal 
splitter or coupler as part of the network would not perform 
as well. 

An experimental measurement of the effectiveness of the 
matching network in Figure 2 would require standard S- 
parameter measurements by a network analyzer at the input 
port of the circulator. Since there is only one source, no 
averaging is needed over phase differences, as described in 
Part I, Section II-C. By the definition of r{uj), and because 
the network is (theoretically) lossless, small values of r(w) 
in 2-4 GHz imply that the power from the source is being 
accepted by the antennas. 

We also note that for two sources (M = 2), we would be 
able to achieve only half the bandwidth for the same —14 
dB threshold. Clearly, the matching network for two sources 
would differ considerably from the one presented in Figure 
2 and could be reciprocal. We do not design such a network 
here. 

B. Two coupled dipoles 

To illustrate the effect of coupling between two parallel 
dipoles, we examine their bandwidth as a function of the 
distance between them. We assume each dipole has the same 
2.4 GHz half-wavelength structure as Figure 1(a). The center 
feeding points of the parallel dipoles are at the same vertical 
level, and the distance between them is d, as measured at any 



Fig. 2. Matching network of the type in Figure 5{b) in Part I is used to 
connect a source (on the left) to two decoupled dipoles (connected on the 
right, typically through baluns that are not shown). Matching networks 1 
and 2 are designed using the Real Frequency Technique [19], followed by 
realizations using Darlington synthesis [20]. 



Fig. 3. Bound 3 for two parallel dipoles versus spacing d, for d between 
0.03A and 1.5A (A =125 mm is the wavelength at 2.4 GHz). Each dipole has 
the structure shown in Figure 1(a), the center feeding points of the dipoles are 
at the same vertical level, and the distance d is maintained along their lengths. 
The horizontal dotted lines indicate Bound 3 for two decoupled dipoles. Also 
shown are the even- and odd-mode beamforming bounds. 


place along their lengths. For each value of d, we apply the 
modeling recipe for the S-matrix detailed in Section II 

in the range 1-5 GHz. Since we wish to compare the bounds 
for coupled dipoles with (10), we model their S-matrices using 
six poles and six zeros, and enforce Sl{oo) = —I. Then the 
resulting model Sl{s) satisfies (1) at Sq = cxd for every d, and 
Bound 3 can be applied to compare with (10). 

We let d range from 0.03A to 1.5A, with step size O.OIA, 
where A = 125 mm is the wavelength at 2.4 GHz. The poles 
and zeros PL,i, are computed from Sl{s), and Bound 3 is 
computed for each d. For M = 1 and M = 2, the bounds are 
shown in Figure 3 by the black and blue curves. The black 
curve values are twice the blue curve. Horizontal dotted lines 
indicate the bounds for decoupled dipoles. 
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Figure 3 shows oscillatory behavior and suggests that cou¬ 
pled dipoles can have large bandwidths at certain distances 
from each other. The maxima appear near 0.24A and are 20% 
larger than their decoupled-dipole counterparts. As d increases, 
the bounds approach the decoupled-dipole limits. On the other 
hand, when d approaches zero, the two dipoles merge into a 
single dipole, and the bounds in Figure 3 approach the bounds 
where M sources drive a single dipole. 

We already know from Part I, Section V-A, that using 
these dipoles in a beamformer conhguration with M = 1 
generally cannot achieve the bound for any d. But the beam- 
former configuration is still worth analyzing briefly. There 
are two beamformers that are of special interest; even-mode 
and odd-mode, corresponding to v{s) = [l/-\/2, l/v/2]^ and 
[ 1 / 72 ,- 1 / 72 ]^ in (25) in Part I. Since the dipoles are 
reciprocal and have a symmetric structure, the resulting Sl{s) 
is symmetric circulant, and both even and odd v{s) are real 
constant unit eigenvectors of Sl{s). Thus Theorem 5 in Part I 
applies. 

Figure 3 shows the results. Both even and odd-mode bounds 
are strictly smaller than the bound for M = 1, as expected. 
Their sum, however, equals the M = 1 bound because 
the incident and reflected waves to and from the loads are 
orthogonal for the two modes. Therefore, the modes can be 
treated as decoupled loads, each with its own equivalent S- 
parameter and bound. 

C. Multiple models for the same loads 

For a set of loads with a measured more than one 

model Sl{s) may be created within a given error tolerance. 
To illustrate this, we take the parallel dipoles in Section III-B 
with measured in 1-5 GHz as examples. For each d, 

we model the S-matrix using the Matrix Fitting Toolbox with 
six poles and six zeros and enforce iS'l(O) = /. The resulting 
Sl{s) then satisfies (1) at Sq = 0, and Bound 1 is applied 
to the coupled dipoles. We contrast the results with Section 
III-B, where Sl{oo) = —I is enforced and Bound 3 applies. 

For M = 1, Figure 4 compares Bound 1 (left y-axis) with 
Bound 3 (right y-axis) for the parallel dipoles when d varies 
from 0.03A to 1.5A. The horizontal dotted line indicates the 
limiting values as the dipoles become decoupled. Both curves 
show similar trends when the dipoles are closely spaced, and 
suggest a bandwidth peak at 0.24A. For d > 0.32A, the curves 
seem out of phase with each other. Further exploration is 
needed to see if there is any physical significance to this. 

D. Rational approximation of time delays 

Any characterization of physically separated antennas 
should account for the time of propagation of the signal be¬ 
tween the antennas. Time delays associated with this distance 
can be modeled as a non-rational (exponential) function of s. 

A simple idealized analytical model of time delay is shown 
in Figure 5, which consist of two RC loads and a coupling 
branch between them. The coupling is capacitive and includes 
a delay component. We examine our ability to match to this 
load in the 100-400 MHz band. As a reference value, without 


X lo‘“ X 10'‘“ 



Fig. 4. For one source and two parallel dipoles with spacing d. Bound 3 
(left jz-axis) is compared with Bound 1 (right j/-axis). Bound 3 is the same 
as the black curve shown in Figure 3. The dotted horizontal line indicates the 
limiting values for decoupled dipoles. 


50pF= 

50 £2 

20 pF= 

0.5 ns ] 

^ 50£2 

50pF= 



Fig. 5. Two capacitively coupled loads. The coupling includes an ideal delay 
component of 0.5 ns. 


the time delay (zero delay), the S-matrix of the loads is real- 
rational and satisfies Sl{oo) = —I, has two poles at —8.00 x 
10® and —4.44 x 10®, and has a zero at 0 with multiplicity 
two. We apply Bound 3 and obtain B = 9.77 x 10®. 

In order to apply the broadband bounds to Figure 5, we 
need to approximate the delay using a rational function. We 
consider two ways of doing this. The first way models the time 
delay term ® using the Fade approximation. We use 

the fraction 

s'^-2xlO^°s^-H.8xlO^°s^-8.4xlO^°s+l.68x10^^ 

s<‘+2x10U>s3 + 1.8x1020s2+8.4x1029s-|-1.68x103S> ' 

This fraction is then inserted in place of the time delay, thus 
making the rational matrix Sl{s) that satisfies Sl{oo) = —I, 
and has poles and zeros shown in Table 11(a). (We omit the 
expression of Sl{s) because of its complexity.) Compared 
with the true S-matrix of the loads, the rational S-matrix has 
an average error of —119 dB in its entries between 100 MHz 
and 400 MHz. We apply Bound 3 and obtain 1.24 x 10®. Using 
(9) with T = 0.2, we obtain SB « 9.64 x 10^, which is 8% 
of the bound. Hence (6) is 
^4x10® 1 

/ log —< 1.34 X 10®. (12) 

/ixio« r'{uj} 

This bound is larger than without the delay. 

A second way to model delay is to avoid the analytical 
Fade step, and directly fit numerical simulations using rational 
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TABLE II 

Values 0¥pL,i and ZL,i for Sl{s) model for the coupled loads 
IN Figure 5. Table (a) shows and ZL,i when the Fade 
approximation is applied, and table (B) shows pL,i and when 
the loads are fitted with rational matrices without the Fade 
approximation. 


M. 


i 

PL.i 

2L.i 

1 

-6.71 X 10« 

0 

2, 3 

(-0.05 ± 3.37i) X 109 

±3.38i X 10® 

4, 5 

(-O.OOit 1.33i) X 10i° 

±1.33) X 10^0 

6 

-4.35 X 10® 

0 

7, 8 

(-0.16± 1.25i) X 10® 

±1.32i X 10® 

9, 10 

(-0.02±6.81i) X 10® 

±6.81i X 10® 


M 


i 

PL.i 

ZL.i 

1 

-6.71 X 10® 

0 

2, 3 

(-0.05 ± 3.37j) X 10® 

±3.387 X 10® 

4, 5 

(-O.OOit l.Olj) X lOio 

±1.01i X lOio 

6 

-4.35 X 10® 

0 

7, 8 

(-0.16± 1.25i) X 10® 

±1.327 X 10® 

9, 10 

(-0.02 ± 6.76i) X 10® 

±6.767 X 10® 


matrices. To illustrate this, we treat the loads in Figure 5 as a 
black box, obtain a numerical frequency response, and use the 
Matrix Fitting Toolbox to fit this response. We use ten poles 
and ten zeros to fit the loads in 100-400 MHz, and enforce 
Sl{oo) = —I- The resulting Sl{s) has poles and zeros shown 
in Table 11(b); the average error between the entries true S- 
matrix of the loads and entries of Sl{s) is —150 dB. Bound 
3 is applied to the fitted rational model, which yields B = 
1.24 X 10®, and 5B « 1.75 x 10^, which is 1% of the bound. 
The total is B + 5B = 1.26 x 10®. This numerically fitted 
model result is similar to the Fade approximation result (12), 
illustrating the robustness of this approach. 

In the next section, we further illustrate the numerical 
method by fitting loads that include sections of realistic 
transmission lines. No analytical modeling is needed. 

E. Four commercial 2.5 GHz antennas 

We consider a system where M = N = A, and the 
loads consist of two pairs of Skycross iMAT-1115 commercial 
antennas designed for 2.5 GHz. The antennas are shown in 
Figure 6(a), numbered 1 through 4 from left to right. These an¬ 
tennas are simulated using Ansys HFSS in 2-4 GHz, and their 
de-embedded S-matrix entries ai'e plotted in Figure 

6(b). The de-embedding removes the effect of the grounded 
coplanar feed lines that drive the antennas in the simulation. 
Each antenna in Figure 6(a) is one-tenth of a wavelength 
away from its neighbor at 2.5 GHz and the antennas therefore 
couple. From Figure 6(b), we see that antennas 1 and 2 are 
decoupled at 2.5 GHz, but not at other frequencies; antennas 
2 and 3 have significant coupling near 2.5 GHz. We wish to 
examine the achievable bandwidth performance of these four 
antennas. 

We apply the modeling method in Section II to the simulated 
Figure 6(b) compares the simulated with 

the S'l(s); the expressions for Sl{s) are not shown. The 
maximum of any element of SlUu)) — in 2-4 GHz 

is -38 dB, and the average over the frequency band is -53 



Fig. 6. (a) Geometries of two pairs of Skycross iMAT-l 115 antennas, (b) The 
simulated and modeled S-matrix elements of the antennas versus frequency 
in 2-4 GHz. The HFSS simulated result S£(ya;) is shown in blue, and 
are compared with S^ljui) in red. Different elements are distinguished by 
different markers in the plot; only some of the elements in the S-matrix are 
shown because of the symmeti'y in the structure of the antennas. 


TABLE III 

Values oe and z^^i eor the S-matrix oe the antennas in 
Figure 6(a). 


i 

PL.i 

ZL.i 

1 

-1.54 X 10® 

1.64 X 10® 

2 

-7.74 X 10® 

3.67 X lO^® 

3, 4 

(-0.43 ± 1.377) X lO^o 

(-0.23 ± 2.237) X lO^® 

5, 6 

(-0.15 ± 1.607) X lO^® 

(-0.22 ± 1.697) X 10^® 

7, 8 

(-0.08 ± 1.697) X lO^o 

(0.04 ± 1.707) X 10^° 

9, 10 

(-0.75 ±2.407) X lO^® 

(2.65 ± 3.467) X 10^° 

11, 12 

(-0.40 ±3.147) X lOi® 

(0.03 ± 1.637) X lO^o 


dB. The pL^i, for the antennas are listed in Table III. It 
is readily checked that det5'L(0) = 1, and (1) is satisfied at 
So = 0. Bound 1 yields 


pUJ2 1 

/ log — -dw < 2.31 X 10“^®, (13) 

icai f’(w) 

where wi = drr x 10® and W 2 = Stt x 10®. 

We apply the steps in Section II-B to estimate SB for r = 
0.2. The integral in (9) with /(w) = 1/w^ yields SB 7.88 x 
10“^^, which is equal to 34% of the bound in (13), and which 
we accept. Therefore, 



—-—-div < B + SB = 3.10 X 10 
r'(w) 


(14) 


is our bound for the four antennas. 

The right-hand side of (14) is compared with the achieved 
bandwidth of two different matching networks. The first net¬ 
work is shown in Figure 7(a) for a single antenna; this network 
is duplicated four times, one for each antenna, and does not 
account for any antenna coupling. The second network, shown 
in Figure 7(b), is a decoupling network at 2.5 GHz designed 
using Method 4 in [21]. The achieved integral for the two 
networks are 6.98 x 10“^® and 2.49 x 10“^^, respectively. 
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Cll 

9.06 nH 

C47 

0.31 pF 

Cl 5 

0.45 pF 

C48 

0.65 pF 

C22 

4.22 nH 

C55 

3.33 nH 

C25 

0.42 pF 

C 56 

0.07 pF 

C 26 

0.54 pF 

C57 

0.10 pF 

C33 

20.46 nH 

C58 

0.004 pF 

C35 

0.05 pF 

C66 

12.02 nH 

C36 

6.44 nH 

C 67 

39.76 nH 

C37 

0.77 pF 

C68 

0.10 pF 

C44 

5.62 nH 

C77 

2.62 nH 

C45 

19.98 nH 

C78 

0.07 pF 

C46 

122.24 nH 

Css 

3.49 nH 


(b) 




Fig. 7. (a) A two-port matching network for a single Skycross iMAT-1115 

antenna (duplicated four times) that does not account for antenna coupling, 
(b) A decoupling network at 2.5 GHz designed using Method 4 in [21], where 
each line represents a capacitor or inductor, and each port is grounded through 
a component not drawn in figure. The capacitance and inductance values are 
listed in the table to the right; ca i = 1,..., 8 are the components connecting 
port i to ground, and Cij i ^ j = 1,..., 8 are the components connecting 
port i and j. (c) The ‘r(cj) for networks (a) and (b) are shown in 2-A GHz. 
Also plotted is r{uj) for (b) when the reactive components have Q factors 10, 
30 and 100. 



Fig. 8. An RF system where M uncorrelated sources drive N loads having S- 
matrix Sl through a passive (M-l-A^)-port matching network S (the complex- 
frequency argument s is omitted). The N loads ai‘e shown in the dashed box 
using their equivalent Darlington representation, which consists of a lossless 
2A^-port network Si, connected with N characteristic impedances Zq. Sc is 
the (M + N) X (M + A^) S-matrix of the concatenated network of S and 
Sb- Slm^ Sq and Sqm Ibe S-matrices as seen at different ports of the 
system. 


standard rational fitting; We employ the Matrix Fitting Toolbox 
to fit the antennas using a rational matrix with 15 poles and 
zeros, and enforce (1) at Sq = 0. We omit the details, but 
the result of Bound 1 is B = 2.74 x 10“^°. From (9), 5B « 
3.61 X 10“^^, and the result is B + 6B = 3.10 x 10“^°, which 
coincides with (14). Hence, the feed lines in this example 
do not affect the bandwidth, likely because they are equal in 
length and have minimal coupling between them. 

We now provide derivations of the bounds used in both Parts 
of this paper. 


Although the decoupling network has a better bandwidth 
performance than the network in Figure 7(a), there still exists 
a significant gap between the achieved integrals and (14), 
indicating that much better bandwidth performance with these 
antennas is still possible. The r'{u;) for both networks is 
plotted in Figure 7(c), where we see narrow bandwidth for 
both networks. Neither network comes close to achieving the 
desired threshold of r = 0.2 (power reflection ratio —14 dB) 
over 2-4- GHz, which should be possible according to (14). 
The design of a network to achieve (14) remains an open 
problem. 

The reactive components in Figure 7(b) are ideal in that 
they have no resistive properties. To illustrate the bandwidth 
lost when the components have finite Q factors, we let the 
reactive components behave as ideal reactive components in 
parallel with resistances. All of the components in Figure 7(b) 
are assumed to have the same Q factors for all frequency, 
where Q is defined as the ratio between the susceptance and 
the conductance of the components. Figure 7(c) shows the 
r'{uj) for different Q factors. 

Finally, we show that the grounded coplanar waveguide feed 
lines used for these antennas have minimal effect on the bound. 
We modify the S-matrix generated by Ansys HFSS by not de¬ 
embedding it from the transmission lines, thus including them 
in the frequency analysis. This S-matrix is then subjected to 


IV. Derivations of Broadband Matching Bounds 

The bounds are corollaries of Theorems 2 and 3, which are 
presented in Section IV-B. The presentation of the theorems 
requires the Darlington representation for multiport loads, 
which is itself presented in Section IV-A. The bounds are 
derived in Section IV-C. 

A. Darlington equivalent network 

We consider the RF system shown in Figure 8, where M 
uncoiTelated sources drive N loads through a passive (M -|- 
V)-port matching network. The network has an (M -|- N) x 
(M -I- N) S-matrix S{s) partitioned as 

M N 

N V^2l(s) S22{s))' 

In Figure 8, we transform the N dissipative real loads Sl{s) 
into an equivalent lossless real 2V-port network Sb{s) ter¬ 
minated by N isolated characteristic impedances Zq. In [20], 
Darlington first verified that such an equivalent transformation 
is possible for V = 1. The extension to V > 1 is shown in 
[23]-[26]. The fact that I—S"[{—s)Sl{s) has full normal rank 
ensures that there are N resistors in the Darlington network 
[25, III.3.1]. 


































































We partition the 2N x 2N S-matrix Sb{s) as 


if So is finite, or 




<5'61i(s) Sbl2{s) 
Sb 2 l(s) Sb 22 {s) 


where Sbij{s) are Nx N submatrices. Then the necessary and 
sufficient condition for Sb{s) being a Darlington equivalent 
network for a real-rational ^^(s) is that Sb{s) is real-rational, 
Hurwitzian, bounded and para-unitary, and Sbii{s) = 5 ^( 5 ). 
We do not need to know the exact form of the rest of Sb{s), 
only its existence is needed. 

In Figure 8 , let S'c(s) be the (M + N) x (M + N) S-matrix 
of the concatenated network of S{s) and S't,(s), partitioned as 


M N 

a (^ M f Scll{s) S'cl 2 (s)\ 

N V^e 2 l(s) SMs))' 

Let S'gm(s) be the N x N S-matrix seen from the output 
ports of Sb{s), and Sa{s) be the N x N S-matrix seen from 
the output ports of S{s). It follows that Sg{s) = 822 ( 3 ), 
Sgm{s) = Sc 22 {s), and 


Sgm{s) = Sb22{s) 

+ Sb 2 l{s)SG{s){I - SL{s)SG{s))-^Sbl 2 is). (15) 


We use the notation px.t and Zx,i, i = 1, 2,... to represent 
the poles and zeros over the WCP of Sx (s), where Sx is 
any of the S-matrices or submatrices shown in Figure 8 . In 
addition, we use the subscript “+” to denote the zeros or poles 
that are in the RHP, and ” to denote those in the LHP. 

In Figure 8 , let 03 ( 3 ) and 63 (s) be the x 1 incident 
and reflected signal to and from the isolated impedances Zq. 
Because Sb{s) is lossless, the amount of power delivered to 
the loads equals the power delivered to the resistive part of the 
Darlington network. Since 63 (jw) = 0, the power delivered at 
juj is ||a 3 (ja;)|p, where asijuj) = S'c 2 i(jw)ai(jw). So the 
total power lost to reflection and dissipation can be written 
as |jai(ju;)|p — ||a 3 (jw)|p, and the power loss ratio of the 
matching network and the loads becomes 

r^{u}) = 1 - ^tr{S'i^i(ja;)S'c2i(jw)}. (16) 

In the rest of the presentation, we use a particular Darlington 
network, which is given in the following lemma in [9]. 

Lemma 1: There exists a Darlington network Sb{s) such 
that 


Zb 22 +,i — —ZL-,i- (17) 

For such Sb{s), the RHP zeros of Sgm{s) are identical to the 
RHP zeros of S'^(—s) — S'g(s)- 


B. Integral log-determinant of SGAiiji^) 

We now present the theorems on the integral of logarithm of 
det S'gm(jw). We assume that I—Sj) (—s)S'l(s) is full normal 
rank. We also assume (1) is satisfied for some Re{so} > 0, 
and choose positive integer m such that 

/ - SI{-s)Sl{s) = - so)"* + - so)"*+' 

(18) 


/ - Sli-s)SL{s) = A^s-^ + + ... (19) 


if So = 00 . In (18) and (19), Am 7 ^ 0 is defined such that 
the entries of / — SJ){—s)Sl{s) have a zero at s = so with 
multiplicity at least m. When so = jojQ, it is shown in Lemma 
4 in Appendix C that m is even; when sq = 00 , m is also 
even since the left-hand side of (19) is an even function. The 
general broadband matching theorems are as follows: 

Theorem 2: Let S'l(s) satisfy (18) for some Re{so} > 0. 
Then for any passive network S{s) such that I—Sl{so)Sg{so) 
is non-singular, we have 


/ Re[(so + (50 +ja;)“^]log|det 5 ’GM(ja;)|fia; 

-'0 


= Trlog 


det Sl{so) 


IIzC'^O + ZL^i) nz('^0 + ZGM-\-,i) 
n*('S0 - ZL,^) n*(so - ZGM+,i) 


( 20 ) 


and 


pOO 

/ [(so-ja;)-('=+i) + (so+ja;)-('=+i)]x 

Jo 

log I det S'GM(iw)Idw = ^ ^ 

I 

~ ~ So) + f y~l(zGM-r,z — Sq) * 

i i 

- '^{-ZGM+,i - So)”*) 1 


( 21 ) 


I 


for k = 1 ,..., m — 1 , where ±so, isg are excluded in pL,i, 
ZL,i in ( 21 ). If I—Sl{so)Sg{so) is singular, then Re{so} = 0 , 
( 20 ) holds, and ( 21 ) holds for fc = 1 ,..., to— 2 ; for k = to— 1 , 
we let So = jwo, and then have 


pOO 

/ [(wo - w)”™ -f (wo + w)“™] log I det SGMijoj)\duj 
Jo 


> 


(- 1 )”" TT 

TO — 1 




- - JC^O)-^”*-'^ + ( - JCCo)-^"*-') 


~ y~l(~^GM-r,z — jzjp) 


( 22 ) 


Theorem 3: Let Sl{s) satisfy (19). Then 

/ OO 

log I det5'GM(jw)|dw 


> 


2k 


_Ept+E<*+ 2 E 


ZGM+,i 


(23) 


for k = 1, 3,..., TO — 1. Equality in (23) holds if A: 7 ^= to — 1 
or / — Sl{s)Sg{s) is non-singular at s = 00 . 

In Theorem 2, although (18) holds, Sg may still be a zero 
of Sl{s) when Re{so} > 0. This happens if and only if 
Sl{s) has a pole at —sq, which cancels the zero at sg in 
(18). Hence, the poles and zeros of S'l(s) at ±so, isj are 
excluded in the sums in (21). The proof of Theorem 2 appears 


-f ... 
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in Appendix D; the proof uses several preliminary lemmas 
which are introduced in Appendices B and C. 

Theorem 3 is adapted from Theorem 1 in [9], which is 
proven with det SlmUoj) in place of det S'gm(jw) in (23). 
Note Slm{s) is the M x M S-matrix seen from the input 
ports of S{s) (see Figure 8). Although it is assumed in [9] 
that M = N, and that the matching network is lossless and 
reciprocal, and the loads are also reciprocal, these assumptions 
can be relaxed. In fact. Theorem 1 in [9] applies without 
change to non-reciprocal networks and loads; the reciprocity 
of the networks and loads is an unnecessary restriction in the 
model. Our version is needed to handle lossy networks, which 
are not handled in [9]. The proof of (23) follows the same 
arguments used in the proof of Theorem 2 and is omitted. 

C. Proof of bounds 

We relate det ScM{ji^) to r(a;) in (16) using the arithmetic- 
geometric mean inequality: 

= ^tr{/- 5'i^i(jw)S'c2i(jw)} 
>detiI-S,2iiju;)Sg^ijio)y/^ 

> det(S'GM(jw)5'GM(jt*^))^^“ = I detS'GM(iw)p/“- 

The first equality holds if and only if the eigenvalues of 
I — S^i{juj)Sc 2 i{j<^) are all equal; this is equivalent to 
- SaSLr^il - SSSL)iI - SgSl)-^S 2 i having 
equal singular values for all juj, which is Condition 2 for 
equality in [1]. The second equality holds if and only if S{s) 
satisfies 5'c2i ( 5 ) 5^21 (-s) + SGM{s)SQf^{-s) = I. Since 
S't,(s) is para-unitary, this is equivalent to ^2 i(s)52^i(-s) + 
Sg{s)S'q{—s) = /, which is Condition 1. 

Taking the logarithm on both sides of the inequality yields 

logr(a;) > (1/M) log | det S'gm(jw)|- (24) 

When So = jwq, we apply (24) to (21) and (22) for fc = 1, and 
omit - since it 

is non-negative. From Lemma 1, the RHP zeros of Sgm{s) 
are identical to the RHP zeros of S'£{—s) — ^g)®); thus 
Condition 4 is necessary and one of the sufficient conditions. 
This finishes the proof of Bound 1. 

When Re{so} > 0, we apply (24) to (20), and then omit 
ndso-tzGM+,.) jj modulus no smaller than one. Note 

lL(S0-ZGM + ,i) 

Re[(so — + (so + in (20) is positive for any 

Re{so} > 0 and uj. The result is Bound 2. 

We combine (24) and (23) for k = 1, and then omit 
ZGM+,i since it is non-negative. The result is Bound 3. 


V. Conclusions and Future Work 

We have presented a bandwidth analysis for multiport 
matching that applies to an arbitrary number of sources 
and coupled loads, using broadband matching bounds on the 
integral of a power loss ratio. Part I presented the definitions 
and bounds and applied them to settings where the loads could 
be described analytically using rational functions. Conditions 
were given for when the bounds could be met with equality. 


Part II focused on realistic loads, including antennas, and the 
rational fitting needed for accurate bandwidth calculations. 
Proofs of all results were also presented. 

Part I demonstrated that the bound scales generally as N/M 
for M sources and N loads. This scaling is not affected 
by coupling as long as it is not “too strong”. Hence, large 
bandwidth is possible even if the loads are coupled or closely- 
spaced. Some realistic examples in Part II showed how large 
bandwidths can be attained in practice. Although we touched 
upon techniques to attain the bounds, the general practical 
design problem remains open. 

For loads, such as antennas, whose frequency responses 
are available numerically through simulation, rational fitting 
of the S-matrix is required to apply the bounds. We showed 
how the rational matrix is then used to find a bandwidth 
bound for the system being fitted. Increasing the polynomial 
orders of the numerators and denominators of the rational 
functions in the S-matrix generally improves the accuracy 
of the system approximation, but can cause the problem of 
over-fitting, which results in a loose bound. Over-fitting is an 
important issue that deserves further attention. 

Our examples included antennas that are matched near their 
resonant frequency. It would be interesting to see how the 
bounds characterize the bandwidth for electrically small an¬ 
tennas, which are being used below their resonant frequencies 
and are generally considered narrowband. The classical Chu 
bound [28] and some recent advances [29]-[31] provide other 
methods to describe the bandwidth for small antennas. One 
future area of study could be to reconcile these methods with 
bandwidth results calculated using the modeling techniques 
and bounds contained herein. 

The communication-theoretic implications of using broad¬ 
band multiport networks in wireless systems need study. We 
have been concerned primarily with the aspects of network 
design that ensure efficient power delivery to loads such as 
an array of multiple coupled antennas. However, the choice of 
network can also affect the radiation efficiency and far-held 
pattern of an antenna array. Hence, a complete system design 
should consider total power transfer from the transmitter 
amplihers to a far-held receiver. We also did not examine 
the implications of using the wideband matching networks for 
antennas that are intended for both transmission and reception. 

Although we have treated the case where uncorrelated 
sources are driving coupled loads, the reverse situation where 
coupled sources are driving decoupled loads needs separate 
analysis. Such a system could arise when the sources are 
closely-spaced receiver antennas connected to isolated low- 
noise amplihers. The noise hgures of the amplihers would 
play a role in the criterion for determining bandwidth. 

Maximizing the data rate attainable with a prescribed set 
of N antennas is a potentially interesting problem. Conven¬ 
tional narrowband Shannon theory says that MIMO, with 
M independent data streams (sources) in a rich scattering 
environment, achieves rates linear in N when M = N [32]. 
But, as we have seen, the bandwidth attainable for M = W 
is only 1/N of that attainable for M = 1. Thus, MIMO 
with N streams achieves 1 /N the bandwidth of MIMO with 
a single stream. Since the transmission data rate is directly 
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proportional to bandwidth, one could therefore conjecture 
that, with N antennas, one stream is as good as N. This 
bandwidth/multiplexing trade-off needs further study. 
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Appendix 

Theorem 1 is proven in Appendix A. We prove Theorem 
2 in Appendix D using several preliminary lemmas that are 
introduced in Appendices B and C. 


A. Proof of Theorem 1 

Let e{(jj) = r^(w) — Then 

/(w) 


/(w) log = 

r{Lj) 


log 1 + 


/(w) 


log 1 + 


l — e{uj) 




g(^) 

r'2(w) 


duj. 


We show that £(w)/(l — r'^(a;)) < p{oj). 
From (5) in Part I, 


du) 


r^(a;) = 1 — 

u{Sg{I - SgSl)-^{I - S^SlW - SgSl)-^S2i} 

M 

and a similar relation between r'‘^(uf) and S'j^{juj) holds. For 
simplicity of the presentation, we omit the argument juj for 
the S-matrices in the remainder of the proof. We perturb 
+ 5SL{ju)\ a first-order expansion yields 

e{(jj) = r^(w) — « 

- - SaSfr^A.il - SgS^)-^S2i}, 

where Ai is the Hermitian matrix 


=5S^iI-S'^SG)-^{Sld -S'^) 

+ iSG-Sf)iI-S'LSG)-^SSL. 

Because Ai is Hermitian, the matrix S 2 i{I—SGS'jf)~^Ai{I— 
SgS'i)~^S 21 is also Hermitian. Some manipulations of the 
trace give us 

£(cc) < ^tr{^2^(/ - SgS'^)-^{I - SgS'^)-^S2i}, 
where CTi.max is the maximum singular value of Ai. Similarly, 




= j^ir{S^,{I-SGS'^)-^{I-S'^S'L){I-SGS'^)-^S2i} 


> 


1 - a'i 




tr{52" (/ - SgS'^)-^ (/ - 


Hence, 


l-r'2(a;) “ I - 

We now seek an upper bound on which depends on 

the maximum singular values of both 6Sl and {Sg — S'j^){I— 
S'j^Sg)~^- The maximum singular value of 6Sl is <Js^max{oj)- 
To obtain the maximum singular value of {Sg — — 

S'i^Sg)~^, we try to prove that 

I + «„ax(w) - - S'^Sg)-^{I - S'^Sg)-^ 

-{I- S'lSg)-^{SE - S'^){Sg - - S'^Sg)-^ 

(25) 

is a positive semidefinite matrix, which means the maximum 
eigenvalue of the positive definite matrix I + ~ 

is larger than or equal 
to the maximum eigenvalue of the positive semidefinite matrix 
(/- S'^Sg)-^ {SB - S'l){Sg -S'1^){I-S'^Sg)-\ Then the 
square root of the maximum eigenvalue of / -f — 

^L%in(‘^))(-^-'5'L>S'G)"^(/-S'j^S'G)"^ is larger than or equal 
to the maximum singular value of {Sg — S'jB){I — S']^Sg)~^- 
To prove that (25) is positive semidefinite, we simplify it 

I + «™ax(^) - - S'^Sg)-^{I - S'^Sg)-^ 

-{I- s'lSg)-^{sB - sf){SG - s'jB){i - S'M-^ 

= {I - S'lSg)-^A2{I - S'lSgTB 

where 

^2 = - <..„(^))/ + {I- S'r.Sf) 

-sB{i-s'iBs'^)SG. 


A 2 is a positive semidefinite matrix because the singu¬ 
lar values of S'^ and Sg are no larger than one, and 
therefore the minimum eigenvalue of the positive definite 
matrix {cr'l^^{uj) - + (/ - S'^S'B^) is larger 

than or equal to the maximum eigenvalue of the positive 
semidefinite matrix Sq{I — S'jBS' j^)Sg- Therefore, we have 
proven that the square root of the maximum eigenvalue of 

I + {^'LM - ^'lM){i - s'^Sg)-^{i - s'^Sg)-^ 

is larger than or equal to the maximum singular value of 
{Sg — S'jB){I — S'j^Sg)~^- The maximum eigenvalue of 

I + «.ax(^) - <^IM){I - SiSGr^{I - SiSG)-^ 


is smaller than or equal to 1 


yH\ 


Therefore, 


the maximum singular value of {Sg — S'jf){I — S'j^Sg) 


-1 


smaller than or equal to W 1 -f 


( 1-0 


xM)" 


-. We conclude 


^l.max < 2(75 ,max (w)l/l -I- 




( 1 - 


L,max 




which yields e(w)/(l — r'^(a;)) < p{oo). This finishes the 
proof. ■ 


M 
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B. Preliminary lemmas on real-rational functions 

Lemma 2: Let /(s) be a real-rational function with /(sq) = 
c 7 ^ 0. Then the series expansion of log/(s) around s = Sq 
can be written as 


log/(s) = logc + ai(s - So) + a 2 (s - Sq)^ + ■ • ■ 

+ — SqY + ..., (26) 


where 




{.Pi -So) ^ - '^{Zi - So)' 


(27) 


for £ = 1,2 ,..where pi, Zi are the poles and zeros of /(s). 

Proof: Since /(so) 7 ^ 0 is finite, log/(so) = logc is 
finite. So the expansion of log/(s) around s = so can be 
obtained using Taylor series expansion. The result is (27). ■ 
Lemma 3: Let /(s) be a real-rational function. For any 
Re{so} > 0 , if /(so) is finite, non-zero, and 1 — /(—s)/(s) 
has a zero at s = sq with multiplicity m, then 


nOO 

/ Re[(so - jw)"^ -b (so + juj)~Y^og\f{juj)\duj 
Jo 

,ni(so -F-t,i)ni(so + ^:+,*) 


= Trlog 


/(so) 


ni(so+p+.i)ni(so-^:+.*) 


(28) 


and 


pOO 

/ [(so - + (so + ju;)-('=+i)]log|/(ju;)|dcc 

^0 


i i 

- ( ~ - So)”'") 

i i 


(29) 


for k = 1 ,..., m — 1 , where pi, Zi are the poles and zeros 
of /(s) in the WCP, and p+^i, 2 ;+^ are the poles and zeros of 
/(s) in the RHP. 

Proof: Since /(sq) 7 ^ 0, we begin by applying Lemma 
2 to write the expansion of log/(s) as (26) and (27). For 
convenience, we define another real-rational function /(s) as 

,i\iis-p+A) ni(s+2+.o 


/(s) = f{s) 


(30) 


ni(s+p-r.*)ni(s-^+.o’ 

where and z+^i are the poles and zeros of /(s) in the RHP. 
Then /(s) has no poles or zeros in the RHP, and satisfies 
f{-s)f{s) =J{-s)f{s) and |/(jw)| = |/(ju;)|. We apply 
Lemma 2 to /(s) and get 

log/(s) = logc-b ai(s - So) + a 2 (s - so)^ + ... 

+ a^(s — Sq)^ -b ..., (31) 



Fig. 9. The contours for the integrals (35), (37), (39) and (41) are shown in 
(a)-(d), respectively. The contours are in the clockwise direction. The sections 
of the contours are labeled Ci, C 2 , • • • and their detailed descriptions ai‘e given 
in the proof of Lemma 3. 


+ j “ '®o) ^ “ '^ 0 ) j ■ (33) 

We can expand log/(s) and log/(s) at s = Sq in similar 
forms as (26) and (31). Since /(s) and /(s) are real-rational, 
the coefficients for the expansion of log/(s) and log/(s) at 
s = Sg are a\ and aj, respectively. 

The next step is to separate our discussions into four cases; 
So = 0, jojo, (To, (To + jojo- Since the cases are similar to one 
another, we elaborate more on the sq = 0 case than the others. 

7) So = 0.' Since 1 — /(—s)/(s) is an even real-rational 
function and has a zero of multiplicity m at sq = 0 , m must 
be an even integer. We take the logarithm of /(—s)/(s) = 
/(— s)/(s) = 1 -b 0(s"*) and use the expansion (31) at sg = 
0. Because /(s) is real-rational, the coefficients in (33) for 
So = 0 are real. We therefore obtain |c| = 1, 02 = 04 = ... = 
dm-2 = 0 and Im{ai} = Imjaa} = ... = Im{am-i} = 0 . 

For So = 0, (28) is trivial since both sides are zero. Because 
02 = 04 = ... = dm -2 = 0, (29) is also trivial for k = 
2,4,..., TO — 2. We show (29) for fc = 1,3,..., to — 1 by 
taking the contour integral of the function 

s-(fc+i)log/(s) (34) 


where 


C = /(so) = /(so) 


ni(so-p-r,»)n»(so + ^-t,0 

ni(so Tp-r,*) n*(so - z+Y 


and 


at = at 


1 

1 


'^{p+.i-so) ^ 


Xl(“7'+7 - «o) 


(32) 



along the closed curve shown in Figure 9(a). This function 
is analytic in the RHP; it is also analytic on the imaginary 
axis except the origin and possible zeros and poles of /(s) on 
the imaginary axis, which we denote as jujt- Therefore, the 
following contour integral is zero: 

f s-('=+i)log/(s)ds = 0, (35) 

J C1+C2+C3+C4 
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where Ci is the line segment between —jR and jR excluding 
[— £,e] and [j{uje — e),j{oJi + &)]; C 2 is the right semicircle 
with radius R centered at the origin; C 3 is the right semicircle 
with radius e centered at the origin; and C4 includes the right 
semicircles with radius e centered at juji. We evaluate the 
integral of (34) as follows: 

[ s-('=+i)log/(s)ds 

Jci 

= (- 1 )“*“, log 

+ OiR-'^), 


where the integral from —R to R excludes the intervals [—e, e] 
and [oji — e,uji + e]. Furthermore, 


/ log f{s)ds= / log f{Re^^)d9 

JC2 - 277/2 


= o 


f ^ogR 

\ R 


o-(fe + l) ■ 


^ 77/2 


S log f{s)ds = 


IC 3 


= JTrdk + 


' —7r/2 


jise^^)-Hogfiee^^)d0 


(-l)'‘2'2arg(c) 

ke^ 




[ log f{s)ds = Y, 

J C4 £ J — 77/2 

X log/(ja;^ + ee^^)jee^^d6 = 0(£log£). 


Combining these path integrals and letting R ^ 00 , e ^ 0, 
we have 


/ OO 

^-(fe+i) log \f(juj)\duj = (-l)^7rafc. 

-00 


Since /(s) is real-rational, \f{-joj)\ = \f*{juj)\ = \f{joj)l 
we have log |/(ja;)| is an even function for k = 

1,3,..., m—1. Using (33), we get (29) for A: = 1,3,..., m—1. 
This finishes the proof for sq = 0. 

2) -So = jwo-' Since 1 — f{—s)f{s) is a real-rational 
function and has a zero of multiplicity m at s = jujo, it 
also has a zero of multiplicity to at s = —jujQ. We take the 
logarithm of /(-s)/(s) = /(-s)/(s) = 1 + 0 ((s - jwo)™) 
and use the expansion (31) at sg = J<^o- We therefore obtain 
|c| = 1 , Im{ai} = Imjag} = ... = Imjafc} = 0 for odd k 
and k < m, and Re{a 2 } = Re{a 4 } = ... = Rejofe} = 0 for 
even k and k < m. 

For So = jwo, (28) is trivial since both sides are zero. We 
show (29) by taking the contour integral of the function 

[{jcoo - s)-('=+i) + (jcco + s)-('=+i)] log /(s) (36) 


along the closed curve shown in Figure 9(b). This function 
is analytic in the RHP; it is also analytic on the imaginary 
axis except ijwo and possible zeros and poles of /(s) on the 
imaginary axis, which we denote as juJi. Therefore, 

f [{jcoo - s)-('=+i) + ijiOo + 

J C1+C2+C3+C4 


X log f{s)ds = 0, (37) 

where Ci is the line segment between —jR and jR excluding 
[j{ujo-e),j{ujo + e)], [j{-ujo - e),j{-ujo + e)] and [ 7 ( 0 ;^- 
£), j(w^-|-£)]; C 2 is the right semicircle with radius R centered 
at the origin; C 3 includes the right semicircles with radius 
£ centered at the juJo and — jwo; and C 4 includes the right 
semicircles with radius £ centered at juJi. 

By evaluating the integral paths in (37), and letting R —> 00 , 
£ —)■ 0 , we obtain 

/ OO 

[(wo - 0 ;)"^''+^^ -f (wo -f log \f{juj)\du} 

-OO 

= (-l)''j'''^^27rafc. 

We have [(wq — -f (wg + log |/(ja;)| is 

an even function, and using (33) we get (29) for k = 
1, 2,..., TO — 1. This finishes the proof for sg = jwg. 

3) So = CTg.' Since /(s) is a real-rational function, 
Im{di} = Im{a 2 } = ... = Im{dm-i} = 0. For sg = (jg, 
we show (28) and (29) by taking the contour integral of the 
function 


[(ag - + (ag + sY'^+^'>] log f{s) (38) 

along the close curve shown in Figure 9(c). This function 
is analytic in the RHP except ag; it is also analytic on the 
imaginary axis except possible zeros and poles of /(s) on 
the imaginary axis, which we denote as jutg. Therefore, the 
following contour integral is zero for k = 0,1,2,... ,m — 1 : 

[ [(ag - + (ag + 

C1+C2+C3+C4 

X log f{s)ds = 0, (39) 


where Ci is the line segment between —jR and jR excluding 
[j{uje — e), j{uje + e)]', C2 is the right semicircle with radius R 
centered at the origin; C 3 is the circle with radius e centered 
at the (Tg; and C4 includes the right semicircles with radius e 
centered at juje. 

To show (28), we evaluate the integral paths in (39) for 
fc = 0, and let R —> 00 , £ —> 0. The result is 

/ OO 

[(CTg - + (CTg -f ja;)"^]log|/(ja;)|dw 

-OO 

= 27rlog |c|. 


We have [(CTg — juj) ^ -I- (CTg + jw) log |/(jw)| is an even 
function, and using (32) we get (28). 

To show (29), we evaluate the integral paths in (39) for 
k = 1, 2,..., TO — 1, and let i? — 00 , £ —>^ 0. The result is 


/ OO 

[(cto - + (CTg -f ja;)"(''+^^]log|/(ja;)|(iw 

-OO 

= (-l)''27rafc. 


We have [(CTg - ju) -f (CTg -f ju) ('=+ 1 )] log |/(jw)| 

is an even function, and using (33) we get (29) for k = 
1, 2,..., TO — 1. This finishes the proof for sg = ctq. 
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4} So = ctq + juJo- For sq = ao + jojq, we show (28) and 
(29) by taking the contour integral of the following functions 

[(so - s)-('=+i) + (so + + (s* - s)-('=+i) 

+ (4 + log /(s) (40a) 

[(so - s)-('=+i) + (so + - (s5 - s)-('=+i) 

-(sS + s)-('=+i)]log/(s) (40b) 


along the close curve shown in Figure 9(d). These functions 
are analytic in the RHP except so and sjj; they are also analytic 
on the imaginary axis except possible zeros and poles of /(s) 
on the imaginary axis, which we denote as juj£. Therefore, the 
following contour integral is zero forfc = 0,l,2,...,7Ti—1: 

[ [(so-s)-('=+i) + (so + s)-('=+i) 

J C1+C2+C3+C4 

+ (s5 - s)-('=+i) + (sS + s)-('=+il] log f{s)ds = 0, (41a) 

f [(so-s)-('=+i) + (so + s)-('=+i) 

■^Ci +C2+C3+C4 

- (s5 - - (sS + s)-('=+il] log f{s)ds = 0, (41b) 


where Ci is the line segment between —jR and jR excluding 
[j{uji — e),j{uj£ + e)]; C 2 is the right semicircle with radius 
R centered at the origin; C 3 includes the circles with radius e 
centered at the so and Sq ; and C4 includes the right semicircles 
with radius e centered at juj£. 

To show (28), we evaluate the integral paths in (41a) for 
fc = 0, and let R —> 00, e —> 0. The result is 

/ OO 

Re[(so - + (so + jw)"^]log|/(ja;)|fiw 

-00 

= 27rlog |c|. 

We have Re[(so — + (so+jw)“^] log |/(ju;)| is an even 

function, and using (32) we get (28). 

To show (29), we first evaluate the integral paths in (41) for 
fc = 1, 2,..., TO — 1, and let i? —)■ 00, e —>■ 0. The result is 

/ OO 

[(so - + (so + log \f{ju;)\du; 

-OO 

= {-l427rdk. 

We have [(sq — + (so + log |/(jw)| is 

an even function, and using (33) we get (29). This finishes the 
proof for So = CTo + Jujq. ■ 

5) Remark on lemma 3: Because the formulas (28), (29) 
subtract or divide poles from zeros in equal quantity, when we 
apply Lemma 3 to /(s) = det ^(s), (28), (29) still hold if the 
poles and zeros of the determinant are replaced by the poles 
and zeros of the matrix ^(s). 


C. Preliminary lemmas on S-matrices 

Lemma 4: If Rejso} > 0 and the Darlington network given 
in Lemma 1 is used, then 

det Sgm{so) = det 5622 ( 50 ) 4 0, (42) 



{PGM,i — So) 


-i 


y^XzGM,i - So) 


-r 


y~l(Pb22,i — So) ^ — ''^{Zb22,i — So) ^ j (43) 

i i / 

for f = 1 ,..., TO — 1. 

If Rejso} = 0, then to is even, (42) holds, and (43) holds 
for f = 1 ,..., TO — 2. 

If Re{so} = 0 and I — Sl{so)Sg{so) is non-singular then 
(43) holds for f = TO — 1. 

If Rejso} = 0 and I — Sl{so)Sg{so) is singular then 

(-l)^x 

y^XPGM,i - So)”*-’””^^ - y^X^GM,i - So)”*-”"”^^ 

i i 

< (-l)^x 

^(P622,i — So)”*-™”^^ — ^(^622,1 — So)”^™”^^ 
i i 

(44) 

Proof: We separate our discussion into two possibilities: 
Rejso} > 0 and Rejso} = 0. 

1) Re{so} > 0.' Because the Darlington network in Lemma 
1 is assumed, 5622 (s) has no zeros at s = so, for otherwise 
5 l(s) would have zeros at s = —so^ ™d therefore have poles 
at s = So in order to satisfy (18). This contradicts 5 l(s) being 
Hurwitzian. Hence, det 5622(so) f 0, and Sl{s) has no zeros 
at s = —So. 

Since Sb{s) is lossless and therefore S'[{—s)Si,{s) = I, we 
get 5612 (s) = -5)(^(-s)5^;^(-s)5622(s). We manipulate 
(15) to get 

5gm(s) = [/ - 562i(s)5g(s)(/ - 5l(s)5g(s))-i 
X 5^^(-s)5,^2i(-s)]5622(s). 




Taking determinant on both sides yields 
det 5 gm(s) = det Sb 22 {s) det[/ - 562 i(s)5g(s) 

X (/- 5g(s)5g(s))-15^^(-s)5,^2i(-s)] 

= det5622(s) det[/- 5 ^i(-s)562i(s)5g(s) 
x(/-5g(s)5g(s))-15Z^(-s)]. 

Since 5 l(s) and 5 g(s) are bounded, 5 g(s)5g(s) is also 
bounded and I — 5g(so)5g(so) is non-singular for Re{so} > 
0 [24, 7.22]. Hence, 5f2i(-s)562i(s) = / - 5 |;(-s)5l(s) = 
0{{s — So)™), and 5 l(—sq) and I — Sl{so)Sg{so) are non¬ 
singular. We then have 

det5GM(s) = det5622(s)[l + 0((s - so)™)]. 

Thus det 5622 (so) = det5GM(so) f 0, and (42) holds for 
Re{so} > 0. 

To show (43), we apply Lemma 2 to det 5622 (s) and 
det5GM(s): 

logdet5GM(s) = ao -I- ai(s - so) + ... 

+ am-i(s-so)™ ^ + ... 
logdet5622(s) = 5o -f 6i(s - So) + ... 

+ bm-l{s — So)™ ^ -f . . . , 

where a£ and b£ have the form (27). Because det5GM(s) = 
det 5622 (s) + 0((s — So)™), or = 6^ for f = 0,1,..., to — 1. 
Writing a£ and b£ in the form of (27) yields (43). 


and 
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2) /?^{5o} = 0; Let sq = ju^o- Since Sb{s) is lossless and 
(18) is satisfied, S'f (jwo)5'L(jwo) = S^^^{ju}Q)Sb 22 {j^o) = 
I. Hence, det S't, 22 (jwo) 7 ^ 0. 

We begin by showing that m is even. We substitute s = 
j(wo±e) into (18); 

■f- 5'r(-j(wo ±£))5 'l(j'(wo ±£)) 

= I - S^UiuJo ± e))SL{j{uJo ± £)) = £l„(±j£)'” + . . . . 

Since 5 'l(s) is bounded, the £lm(±j£)’” is positive semidefi- 
nite. With Am ^ 0, it is possible only when m is even. 

If I— SL{joJo)Sc{joJo) is non-singular, we follow a method 
similar to the Re{so} > 0 case to get det S'gm(s) = 
detS'h 22 (s) + 0{{s — jwo)™). Thus (42) and (43) hold for 

£ = 1,2,..., TO — 1. 

If I — S'l(jwo)<S'g(jwo) is singular, (44) can be proven in 
a manner similar to the proof of Lemma 5 in [9]. These steps 
are omitted. This finishes the proof of Lemma 4. ■ 


D. Proof of Theorem 2 

We use the Darlington network in Lemma 1. Then Lem¬ 
ma 4 gives detS'h 22 (so) = det5GM(so) f 0. Because 
Sb{s) is para-unitary, (18) implies det(5'^2(~'5)'5'&22('S)) = 
det(5'f (-s)5'l(s)) = 1-I-0((s - so)™). Hence, detS'h 22 (s) 
satisfies the conditions of Lemma 3. 

From Lemma 4, det S'gm(s) = det 5'622(s)-l-0((s —sq)™) 
for Re{so} > 0, and detS'GM(s) = detS' 622 (s) -f 0((s — 
for Re{so} = 0. When Re{so} = 0, we let sq = jujo 
and consider s = j (wq ± s) for £ > 0: 

det[S'GM(-j'(t*^o ± e))5'GM(j(t*^o ± e))] 

= det[5'eM(j'(‘*^o ± e))5'GM(j'(‘*^o ± £))] 

= 1 + 5^_i(±j£)™-i+0(£™)<L 

The inequality is because S'gm(s) is bounded. Since to — 1 is 
odd, bm-i = 0. Hence, det(S'Q^(— s)S'gm(s)) = l-l-0((s — 
So)"*), and det5'GM(s) satisfies the conditions of Lemma 3. 

Unfortunately, Lemma 3 does not apply to det Sl{s) since 
there are possible zeros of S'l(s) at s = sq when Re{so} > 0. 
From (18), if Sl{s) has zeros at sq, it also has zeros at Sq 
and poles at — sq,— sj; the multiplicities of these poles or 
zeros are equal. We therefore construct a function detS'L(s) 
by removing the poles at — sq,— sj and zeros at so,Sq from 
det S'l(s), such that |detS'L(jw)| = | det 5'L(jw)|. This 
function satisfies det S'l(— s)det S'l(s) = 1 + 0((s — so)"*) 
and det Sl(so) f 0, thus the conditions of Lemma 3. 

Since Sb{s) is para-unitary, | det S'L(jw)| = | det S' 622 (jL;)|- 
It follows that 

I detS'L(ja;)| = |detS'L(jw)| = | det S' 622 (jw)|- (46) 


We first apply (28): 

POO 

/ Re[(so - jw)“^ + (so + jw)“^] log| detS'GM(jw)|(iw 
Jo 

= TTlog detS'GM(so) • ^4^7-^ 

1L(S0 - ZGM+,r) 

pOO 

/ Re[(so - jw)“^ -f (so -f jw )“4 log| detS'b22(jw)Mw 
Jo 


(47a) 


= Trlog 


detS'b22(so) 


ni('^0 + Zb22+,i) 
ni('S0 - Zb22+,i) 


(47b) 



juj) ^-f(so-fja;) 4 log |detS'L(jw)|da; 


= Trlog 


det S'l(so) 


ni('So + ZL+,i) 

ni(so - ZL+,^) 


(47c) 


Note So and Sq are excluded in ZL+,i in (47c). From the 
definition of det S'/,(s), we can rewrite (47c) as 



Re[(so-jw) ^ + {so+juj) 4 log |detS'L(ja;)|dw 


= Trlog 


det S'l (so) 


ni(^o+ ^L+,i) 

n*(so - ZL+,i) \ ’ 


(48) 


with So and Sg included in ZL+,i- Because of (46), the 
integral in (47b) is the same if we replace | det Sb 22 {joj)\ with 
I det S'l(j w)|- Hence, the right-hand sides of (47b) and (48) are 
equal. We now apply (42) in Lemma 4 to the right-hand sides 
of (47a) and (47b), and (17) in Lemma 1 to the right-hand 
side of (47b). The result is (20). 

We then apply (29): 


pOO 

/ [(so-jcu)-('=+i) + (so+ja;)-l'=+il] 

^0 

X log| detS'GM(jw)|dw 

= ^ -[ y ^( PGM ,» - So )“^ - '^{zgm,i - So )“^ 

i i 

+ (y^,{ZGM+,i - So)~^ - y^f-ZGM+,i - So)“^) 

i i 

(49a) 


POO 

/ [(so-JCu)-('=+l) + (so+jaT)-l'=+ll] 

7o 

X log I det Sb 22 {.joj)\dijj 

= ^ - |^^(Pb22,i — So)“^ — ^(2&22,i — Sq)~^ 

I i 

+ — So) * — ^( — 2:622+,i — So) 

i i 

pOO 

/ [(so-JCu)-('*+l)+(so+jaT)-l'*+ll] 

^0 

X log \det SL{juj)\duj 

= - so)-'= 

i i 

+ {J2(zl+,^ - so)-" - > 

i i 


(49b) 


(49c) 


where fc = 1, 2,..., to — 1. Note —sg and —Sg are excluded 
in PL,i, and sq and sj are excluded in z^^i and z^+^i in (49c). 
Because of (46), the integral in (49b) is the same if we replace 
I det5&22(jw)| with |det S'l( jw) I - Hence, the right-hand sides 
of (49b) and (49c) are equal. 

When k ox 1 — Sl(so)Sg(so) is non-singular, we 

apply (43) in Lemma 4 to the right-hand sides of (49a) and 
(49b), and (17) in Lemma 1 to the right-hand side of (49b). 
The result is (21). 
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When k = m — 1 and I — S'l(so)<S'g(so) is singular, we 
apply (44) instead of (43) to the right-hand sides of (49a) and 
(49b). The result is (22). ■ 
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